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Abstract 

In this paper, we examine the question of the boundary controllability of the one-dimensional non- 
isentropic Euler equation for compressible polytropic gas, in the context of weak entropy solutions. 
We consider the system in Eulerian coordinates and the one in Lagrangian coordinates. We obtain 
for both systems a result of controllability toward constant states (with the limitation 7 < | on the 
adiabatic constant for the Lagrangian system) . The solutions that we obtain remain of small total 
variation in space if the initial condition is itself of sufficiently small total variation. 
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1 Introduction 

1.1 General presentation 

This paper examines the question of boundary controllability of the non-isentropic Euler equation for 

polytropic compressible fluids in one space dimension, in both Eulerian and Lagrangian forms. The two 

systems under view are the following hyperbolic 3x3 systems of conservation laws, which in our problem 

are considered in a space interval (0,L). First, the usual form of the system, in Eulerian coordinates, is 

as follows: 

d t p + d x (pv) =0, 

d t (pv) + d x (pv 2 + P) = 0, 



t-i„j , D \ , o rr-i 3 



d t ( J -^-p« 2 + Pj+d x [ ^^Pv A + iPv ) = 0. 



(1.1) 



In this system, p = p(t, x) > describes the local density of the fluid at time t € (0, T) and position 
x <G (0, Li), v is the local velocity of the fluid, P > is the pressure. Here 7 > 1 is the adiabatic 
constant. These three equations describe respectively the conservation of mass, momentum and energy. 
In particular the specific total energy E of the system is described as 

E=^v 2 + e, 

the internal energy e being connected to the pressure P by 

P 



^ (7-iy 

In Lagrangian coordinates, the system is written as follows: 




(1.2) 



(1.3) 



Here r := 1/p is the specific volume and e is consequently written as e = -^zr- This system is obtained 
from (1.1) through the change of variable 

y= [ p(t,s)ds, (1.4) 

Jx(t) 

where x(t) is a time-dependent path satisfying 

x'(t) = v(t,x(t)). 

Regular solutions of (1.1) and of (1.3) are equivalent through the above change of coordinates, but this 
turns out to be true even in the case of weak (entropy) solutions (see Wagner [41]) that are under view 
in this paper. However the controllability problems described below are different, since they occur in the 
fixed space domain (0, L), with boundary controls. This domain is not invariant through the change of 
variables (1.4). 

Of particular importance in the study of compressible fluids is the physical entropy function. Setting 
without loss of generality the usual coefficient c v to 1, this function of state reads: 

Regular solutions of (1.1) and (1.3) then satisfy respectively the systems 

d t p + d x (pv) = 0, 

d t v + vd x v + (d x P)/p = 0, (1.6) 

t S + vd x S = 0, 



and 



d t r 


-8yV = 


= 0. 


d t v 


+ d y P-- 


= 


d t S 


= 0. 





(1.7) 



However, in the context of weak entropy solutions, Systems (1.6) and (1.7) are no longer equivalent to 
(1.1) and (1.3). 

Before describing the problem in details, let us specify the type of solutions that we consider. As is 
well-known, these two systems belong to the class of nonlinear hyperbolic systems of conservation laws 

U t + f(U) x = 0, /:SlcMMr, (1.8) 

satisfying the (strict) hyperbolicity condition that at each point df has n distinct real eigenvalues 
Ai,...,A„. These scalar functions are the characteristic speeds at which the system propagates. It 
is classical that in such systems, singularities may appear in finite time even if the initial condition is 
smooth. Hence it is natural from both mathematical and physical viewpoints to consider weak solutions, 
in which discontinuities such as shock fronts may appear. But since uniqueness is in general lost at 
this level of regularity, one has to consider solutions that satisfy additional entropy conditions aimed 
at singling out the physically relevant solution. This paper focuses on entropy solutions with bounded 
variation. The initial state will be supposed to have small total variation as in the framework of Glimm 
[24]. 

We investigate these two systems from the point of view of control theory, and more precisely we 
consider the issue of controllability through boundary controls. This problem is to determine, given an 
initial state of the system uq = (po,vo,Po) or uq — (jo,vo,Po), which final states u\ can be reached at 
some time T > by choosing relevant boundary conditions at x = and x — L (given a notion of such 
boundary conditions). We emphasize that in our problem, boundary conditions on both sides of the 
domain can be prescribed. However the question of determining exactly the set of reachable states seems 
very difficult, since the nature of the system suggests that one has to require additional conditions on u\ 
for it to be reachable (this is in particular connected to an effect of these systems known as the decay of 
positive waves, see in particular Bressan's monograph [7].) Here we will concentrate on the question of 



controllability to constant states. In other words, we aim at proving that given an initial data uq in some 
functional class, it is possible to find a solution bringing the state to a constant. Moreover, we would 
like to focus on the property that the solution should remain of small total variation whenever the initial 
data is of small variation. 

1.2 Mathematical framework 

As mentioned before, we consider in this paper weak entropy solutions, which may present discontinuities, 
in particular shock waves. Let us describe exactly this type of solutions by recalling the basic definitions. 

It will be useful to work with both conservative variables U = (p,pv,pE) and U = (t,v,E) (respec- 
tively for Systems (1.1) and (1.3)) and primitive variables u — (p,v,P) and u — (t,v,P). 

The solutions that we consider will be of bounded variation in space uniformly in time, that is in the 
space L°°(0,T; BV(0, L)) and will not meet the vacuum in the sense that p will be strictly separated 
from (and r bounded). The regularity will be automatically valid for both conservative and primitive 
variables since BV is an algebra. Using the equation leads to the fact that these solutions have a time- 
regularity of class £ip(0, T; L 1 (0, L)). We will denote $7 the domain where the solutions live. It is given 
by {(p,v,P) J p > and P > 0} for System (1.1), and by {{t,v,P) / r > and P > 0} for System (1.3). 
With a slight abuse of notations, we will write U £ ft for the conservative variables as well. 

Now we can consider weak solutions of (1.8) in the sense of distributions; but as mentioned before 
we have to add entropy conditions to the solution in order to retrieve the correct solution. First, recall 
that an entropy/entropy flux couple for a hyperbolic system of conservation laws (1.8) is defined as a 
couple of regular functions (n, q) : il — > R satisfying: 

yUeQ, Dr}(U) ■ Df(U) = Dq(U). (1.9) 

Of course (77, q) — (±Id, ±f) are entropy/entropy flux couples. Then we have the following definition: 

Definition 1. A function U £ L°°(0,T; BV(0, L)) n £ip(0,T; L 1 (0, L)) is called an entropy solution of 
(1.8) when, for any entropy/entropy flux couple {rj,q), with r\ convex, one has in the sense of measures 

V(U) t + q(U) x <0, (1.10) 

that is, for all (p £ X>((0, T) x (0, L)) with ip>0, 

(rj(U(t,x))<p t (t,x)+q(U(t,x))tp x (t,x)) dxdt > 0. (1.11) 

(0,T)x(0,L) 

Now we notice that in Definition 1, we do not mention boundary conditions, which are however of 
primary importance since they compose the control in our problem. Boundary conditions for hyperbolic 
systems of conservation laws are a tedious question, especially when considering entropy solutions. A 
precise meaning of such boundary conditions can be given, see in particular Dubois-LeFloch [20] , Sable- 
Tougeron [36] and Amadori [1]. However in order to avoid this issue, we will rephrase the problem into 
an equivalent one which does not use boundary conditions explicitly. We fix an initial condition as above, 
and consider (1.1) and (1.3) as under-determined systems (without boundary conditions). The question 
is to determine for which states u± there exists a solution u = (p, v, P) or u — (r, v, P) in (0, T) x (0, L), 
with initial state Uq and with U\ as final state at time T . The corresponding boundary values can then 
be retrieved by taking the corresponding traces at x = and x = L. At the level of regularity that we 
consider, this is not problematic. 

1.3 Results 

The two results that we establish in this paper are the following. We begin with the result concerning 
the system in Eulerian coordinates. 

Theorem 1. Let uq :— (p q ,vq,Po) £ M 3 with p ,Po > 0. Let 77 > 0. There exists e > such that for 
any uq = (po,vo,Po) £ BV(0,L;M. 3 ) such that 

\\u - «b|U-(o,z) + TV(u ) < e, (1.12) 



for any u\ — (p 1 , Vi, Pi) with p x , P\ > 0, there exist T > and a weak entropy solution of System (1.1) 
u£ J L oo (0,T; J BV A (0,L))n£ip([0,T]; J L 1 (0 I i )) su ch that 

u\t=o — u o an d uu = t =Ui, (1-13) 

and 

TV(u(t,-))< V , Vte(0,T). (1.14) 

Our second result concerns the system with Lagrangian coordinates. This result is different from at 
least two viewpoints: the range of admissible 7, and the role played by the physical entropy. 

Theorem 2. Suppose that 7 £ (l, |j. Let r\ > 0. Let uq :— (to,Vo,Pq) £ R 3 with tq,Pq > and let 
U\ = {ji ,Vi,Pi) with T\,P\ > 0, such that 

S(u!) > S(u ). (1.15) 

There exist e > and T > such that for any uq = (tq, vq, Pq) £ BV(0, L; R 3 ) such that 

IK - «o|U«(o^) + TV(u Q ) < e, (1.16) 

there is a weak entropy solution u £ L°°(0,T; BV(0,L)) r\Cip([0,T]; L 1 ^, L)) of System (1.3) such that 

«|t=o = u a and u\ t=T = ui 1 (1-17) 

and 

TV(u(t r ))<r,, VtG(0,T). (1.18) 

Remark 1.1. Since —S is a convex entropy for (1.3) in the sense of (1.9), associated with the entropy 
flux q — 0, the condition (1.15) is necessary (or at least, the non-strict inequality is). 

Remark 1.2. We conjecture the result to be false for 7 > | in the same spirit as in Bressan and 
Coclite's paper [8]. See Subsection 1.4 below for a brief description of the result of [8J. 

Remark 1.3. Even in the case where 7 £ (l, |), the controllability of System (1.3) is surprising, due to 
the fact that the second characteristic family of the Lagrangian system has constant characteristic speed 
0. Of course, this is the worst case scenario for boundary controllability, since this means that there is no 
propagation from the boundary to the interior of the domain. Hence one has to rely on the interactions 
of the other characteristic families to act indirectly on the second one. Note that in the context of regular 
solutions of class C , the equivalent result is false, since one cannot modify the physical entropy: see 
(1.7). It is the only example that we know, where there exists a result of boundary controllability in the 
context of entropy solutions, while the equivalent fails in the C 1 framework. 

1.4 Previous studies 

Let us say a few words about previous studies on connected subjects. Questions of boundary control- 
lability of one-dimensional hyperbolic systems of conservation laws have been studied in two different 
frameworks, which give rather different results. 

The first one consists in considering classical solutions of these systems, by which we mean of class 
C 1 ([0, T] x [0, L]). Since such systems develop in general singularities in finite time, the solutions which are 
considered are in general small perturbations in C 1 of a constant state, which ensures a sufficient lifetime 
of the solution for the controllability property to hold. Results of controllability for one-dimensional 
systems of conservation laws and more generally quasilinear hyperbolic systems in this framework of 
classical solutions can be traced back to the the pioneering work of China [12]. Many results of very 
general nature have been obtained in this framework since, see in particular Li and Rao [32] for an 
important work on this problem and the more recent monograph by Li [31]. This framework allows 
to work with very general hyperbolic systems (including those in non-conservative form), the main 
condition being that the characteristic speeds are strictly separated from zero, see again [32], [31] and 
references therein. A result which considers the case of a possibly vanishing (but not identically vanishing) 



characteristic speed is due to Coron, Wang and the author [16]; as we will see, it can be applied to (1.1), 
but not to (1.3); and it considers regular solutions for which the theory is rather different from the one 
considered here. 

The second framework in which the boundary controllability of one-dimensional hyperbolic systems 
of conservation laws has been studied is the one of entropy solutions. One has to underline that the 
situation is very different in both contexts, and not a mere difference of regularity. One of the reasons 
for that is that systems (1.8) cease to be reversible in the context of entropy solutions. The reversibility 
or the irreversibility of a system is of central importance in controllability problems. 

Concerning weak entropy solutions, the study of controllability problems for conservation laws has 
been initiated by Ancona and Marson [3], in the case of scalar (n = 1) convex conservation laws. 
Then Horsin [27] obtained further results on Burgers' equation, by using the return method, which was 
introduced by Coron in [13] (see also Coron's book [14]) and which is also an important inspiration here. 
Another result in the scalar case was recently obtained by Perrollaz [35] when an additional control 
appears in the right hand side. 

In the case of systems (n > 2), controllability issues has been first studied by Bressan and Coclite 
[8]. For general strictly hyperbolic systems of conservation laws with genuinely nonlinear or linearly 
degenerate characteristic fields and characteristic speeds strictly separated from zero, it is shown that 
one can drive a small BV state to a constant state, asymptotically in time, by an open-loop control. 
Another result in [8] , especially important for our study, is a negative controllability result in finite time. 
This result concerns a class of 2 x 2 systems containing the system below (which was introduced by Di 
Perna [19]), and which is close to isentropic dynamics: 




■T" 1 ] =0. 



(1.19) 

The authors prove that there are initial conditions, with arbitrarily small total variation, and for which 
no entropy solution remaining of small total variation for all t, reaches a constant state. An important 
property of System (1.19) to establish this result, is that the interaction of two shocks associated to 
a characteristic family generates a shock in the other family. In particular this allows to prove that, 
starting from some initial data having a dense distribution of shocks, this density propagates over time 
provided that the total variation of the solution remains small. Consequently, the system cannot reach a 
constant state in this case. However this property is not shared by the actual isentropic Euler equation, 
and this was used by the author in [21] in order to establish a result on the controllability of this 2x2 
system. The present paper can be seen as a sequel to [21]. Note that this property of two shocks of a 
family generating a shock in the other family is true for the first and third fields of (1.1) and (1.3) when 
7 > | (at least for weak shocks), but when 7 < | such an interaction generates a rarefaction wave in 
the other family (see in particular Chen, Endres and Jenssen [10]), a fact which is crucial in the proof of 
Theorem 2. The behaviour for 7 > | explains our conjecture of Remark 1.2. We were not able to prove 
estimates of decay of positive waves as precise as in [8] for 3x3 systems; hence a generalization of [8] to 
system (1.3) seems difficult for the moment. 

Other important results in the field are due to Ancona and Coclite [2] , in which they investigate the 
controllability properties for the Temple class systems and to Ancona and Marson [4], in which they 
consider the time asymptotic problem, controlled from only one side of the interval. 

1.5 Short description of the approaches 

The main part of the proof consists in proving the following weaker statements. 

Theorem 3. Let uq : = {Pq^o, Po) & H^ 3 with p , Pq > 0. Let rj > 0. There exist e > and T > such 
that for any uq = (po,Vo,Po) £ BV{0,L) satisfying (1.12), there exist a state u\ with p l7 Pi > and 
a weak entropy solution u G L°°(0,T; BV(0, L)) n Cip([0,T]; L x (0, L)) of System (1.1) satisfying (1.13) 
and (1.14). 

Theorem 4. Suppose that 7 <E (l, f ) . Let n > and S > 0. Let u := (t ,vo, Pq) £ K 3 with Tq,Po > 0. 
There exist e > and T > such that for any uq — (to, Vq, Pq) £ BV(0, L) satisfying (1.16), there exists 



a state u\ £ K 3 with T\,P\ > satisfying 

|«i-«o|<tf, (1.20) 

and a weak entropy solution u £ L°°(0,T; BV(0, L)) n Cip([0,T]; L 1 (0, L)) of System (1.3) satisfying 
(1.17) and (1.18). 

When one has succeeded in reaching one constant state, reaching any constant by remaining of 
small total variation is simple, especially in the case of System (1.1), where this can be seen as an 
immediate consequence of the results [32] and [16] concerning the controllability of hyperbolic systems of 
conservation laws in the framework of regular solutions. System (1.3) having a characteristic field with 
constant zero velocity does not enter this framework though, and the proof needs an additional argument; 
in particular this is where (1.15) intervenes. Precisely, we will show the following two statements. 

Proposition 1.1. Given u a and Ub two constant states in Q and n > 0, there exist T > and u £ 
C ([0, T] X [0,-L]) a regular solution of (1.1) such that u(0, •) = u a , u(T, ■) = uj, in [0,L] and 

||w||c°([0,T]iC 1 ([0,£D) < V- (1-21) 

Proposition 1.2. Given u a and Ub two constant states in fl such that 

S(u b ) > S(u a ), 

and given r\ > 0, there exist T > and u £ C^QOjT] x [0, L]) a regular solution of (1.3) such that 
u(0, •) = u a , u{T, ■) — Ub in [0, L], and 

IMIc°([0,T];Ci([0,L])) < V- (1-22) 

Consequently the main objective of this paper will be to prove Theorems 3 and 4. In both cases, we 
use an idea that was already present in [21], that is to use strong discontinuities. By strong, we mean 
discontinuities that are not intended to be of small amplitude, or to be more accurate that are not seen 
as small. This may seem strange to introduce such material in view of (1.14) and (1.18) in Theorems 1 
and 2. In fact, we will use discontinuities that we will consider large during the main part of the proof; 
our analysis relies on interaction estimates due to Schochet [38]. Only in a final step, we will explain 
why these discontinuities can be taken not so large after all. 

In the case of Theorem 3, the construction relies on a contact discontinuity of the second characteristic 
family, which crosses the domain. Then we use additional waves and cancellation effects to kill the waves 
inside the domain along this strong discontinuity, so that in the end the state in the domain is constant. 

In the case of Theorem 4, the construction relies on two shocks of the first and third characteristic 
families, which cross the domain one after another. In the case of System (1.3), the second family cannot 
be used, having identically zero characteristic speed. The first shock is used to filter some waves, so that 
along the second one one can get rid of the remaining waves, still relying on cancellation effects. 

The method that we employ to construct solutions is an adaptation of the wave front tracking al- 
gorithm, inspired in particular by Bressan's version of the method [6]. It should be underlined that 
there arc other methods to establish the existence of entropy solutions of conservation laws, in particular 
Glimm's random choice method [24] and the vanishing viscosity method, see the very general result of 
Bianchini and Bressan [5] . But we have no idea how to use these approaches in the context of control- 
lability problems for conservation laws. The random choice method can be seen however as a method 
to discretize the control in some cases where the limit system is known to be controllable, see Coron, 
Ervedoza and the author [15]. In the same spirit, the question of being able to pass to the vanishing 
viscosity limit in controllability problems for conservation laws is an active research field, limited for 
the moment to cases where the limit equation is known from the beginning to be controllable, see in 
particular [17, 25, 23, 22, 29, 33]. 

1.6 Structure of the paper 

The paper is organized as follows. In Section 2, we recall some basic tools of the theory of one-dimensional 
hyperbolic systems of conservation laws, and introduce some objects which are needed in the construction. 



In Section 3, we introduce some other objects which are more specific to Systems (1.1) and (1.3). In 
Section 4, we describe the construction for System (1.1). In Section 5, we describe the construction for 
System (1.3). It should be underlined that the construction in the Lagrangian case is also valid for the 
Eulerian case when 7 < I • In Section 6 we prove the convergence of the front-tracking approximations 
constructed in Sections 4 and 5 and conclude the proofs of results presented above. Finally in Section 7 
we make some remarks on the size of the solution and on the time of controllability. 

2 Some tools for systems of conservation laws 

In this section, we recall and introduce some general material for hyperbolic systems of conservation laws 
which is not specific to Systems (1.1) and (1.3). We assume that the reader is familiar with the basic 
theory of one-dimensional systems of conservation laws; we refer to Bressan [7], Dafermos [18], Lax [28], 
LeFloch [30], Serre [39] or Smoller [40]. 

2.1 Notations 

It is be useful to put systems (1.1) and (1.3) in the following form rather than in the form (1.8): 

V(«)t + /(«)* = 0, (2.1) 

where at each point u in the state domain ft 

the matrix — — is invertible. (2-2) 

ou 

This allows to work with primitive variables and to apply the results of Schochet [38]. These systems 
are strictly hyperbolic away from vacuum, that is, 

— — ) — — has n distinct real eigenvalues \\{u) < ■ ■ ■ < A„(u), (2.3) 

ou J ou 

which are the characteristic speeds of the system. To each i = 1, . . . , n is associated the right eigenvector 
Ti, determined up to a multiplicative constant; then we define the eigenvectors in terms of the ip variables: 

and then the corresponding families of left eigenvectors ti and Li which satisfy 

t i -r j =L i -R j = 5 i j. (2.5) 

The systems under view satisfy the property that 

each field is either genuinely nonlinear (GNL) or linearly degenerate (LD), (2-6) 

which corresponds respectively to 

r,; • VA; ^ and to r, • VA^ = in f2, 

In the former case, we will systematically normalize the eigenvectors Vi in order for 

r, ■ VXi = 1 in ft, (2.7) 

to be satisfied. In the latter, we will moreover suppose that 

in the coordinates given by u, the vector field r, is constant with |r;| = 1. (2-8) 

We denote IZi = TZi(ai,u-) the rarefaction curves, that is, the orbits of the vector fields rj. The part 
corresponding to Oi > is composed of points u + which can be connected to u_ from left to right by a 
rarefaction wave. We will refer to couples (u_,u+) with u + = lZi(o~i,U-), o~i < as compression waves. 

8 



We denote Si — <Sj(<Xj, w_) the i-th branch of the Hugoniot locus, which is the set of solutions u + 6 fl of 
the Rankinc-Hugoniot equations: 

f(u+)-f(u-) = a(<p{u+)-<p(u-)), seR. (2.9) 

As usual, given a discontinuity between two states w_ and u+, we write [g] for g(u+) — g{u-) and the 
shock speed is denoted as s — s(ii-,u+). We recall that on Si, one has (see e.g. [18, (8.1.9)]) 

s = X - (Ai(«_) + Xi(u+)) + 0{\u+ - u_| 2 ). (2.10) 

The curve Si is parameterized in order that admissible shocks correspond to negative values of the 
parameter ctj. All along this half curve, these shocks satisfy Lax's inequality 

A,(u + ) < s < A<(«_), (2.11) 

and the discontinuity (u_, u + ) traveling at shock speed satisfies (1.10) (see also (3.25) and Paragraph 6.2.3 
below.) 

We recall that for LD fields, the curves IZi and Si coincide and correspond to states connected to u_ 
via a contact discontinuity (whatever the sign of the parameter) . 

We denote by Tj = Tj(<7,, u) the wave curve associated to the i-th characteristic field. We recall that 
for GNL fields, it is composed of the curves Hi for Oi > and Si for o~i < 0. For LD fields, it is composed 
of the coinciding curves TZi and Si . 

Let us now be more specific about the parameterization of the curves TZi, Si and T t . For the linearly 
degenerate fields, the three curves coincide and one sets 

Ti{o- i ,u)=u + a l r i . (2-12) 

For the genuinely nonlinear fields the curves a parameterized so that (for instance in the case of Tj): 

\{Ti{cTi,u)) - \i(u) = o-i. (2.13) 

This parameterization, with the normalization (2.7), ensures that the wave curves T, are of class C 2 ' 1 

and satisfy 

dTi d 2 T t 

■—^(0,u)=ri(u) and — -^ 

This is a standard computation, see e.g. [7, Section 5.2]. Another important consequence of this choice 
of parameterization is that 



^(0,u)=ri(u) and -^(0, «) = (n ■ V)r,-(«). (2.14) 



u = TZi (—o~,-) oTZi(a, -)u, u = Si(—a, ■) o Si{a, -)u and u = T,(— a, •) o T,(cr, -)u. (2-15) 

We will denote by a = (o^, . . . ,o~ n ) the wave vector of a complete Riemann problem and write 

T(a, ■) := T n (a n , •) o • • • o Ti(cti, •)■ 

It will be useful to use the notation Tj for the wave curve from the right associated to the i-th charac- 
teristic field: 

T;(<7, •) :=Tj(-<7, •), that is, m = T,((Tj,u r ) <=>• u r = T l (a i ,ui), (2-16) 

and 

T(cr, •) := Ti(cri,-) o •• • o T„(cr n , •)■ 

Note that 

T(a,T(a,u))=u. (2.17) 

Solving the Riemann problem consists in two parts. First, given w_ and u+ in CI, one finds a vector <r 
such that 

u + = T{o~, u_). 

This is possible at least when u_ and u+ are close enough (at a distance one from another which is 
uniform as U- lies in a compact of ft) and in that case we denote 

a = S(u_,U-)-). 



Then in a second time, one constructs a self-similar function u(t,x) = U{x/t) as follows. One sets 

Uq = U-, Ui = Ti(<7i,-) o •• • oTi(cti,-)u , 

and determines U by: 

• U(x/t) = U{ for x/t e [Xi(ui),X i+ i(ui)], 

• for x/t £ [Ai(ui-i), Xi(ui)]: when the z-th characteristic field is LD or is GNL and cr, < one writes 
the contact discontinuity/shock: 

x x 

U(x/t)=Ui-i for — < s(ui-i,Ui) and Ui for — > s(ui—i,Ui), (2-18) 

and when the «-th characteristic field is GNL and <jj > one writes the rarefaction wave: 

U(x/t)=K i (a,u i - 1 ) for |=Ai(^(er,«i_i)), cr G [0,^]. (2.19) 

We also recall that a Majda-stable shock [34] is a solution (u_,w+) of the Rankine Hugoniot equations 
(2.9) satisfying moreover that 

s is not an eigenvalue of ( — — ) — — (u±), (2.20) 

\ouJ ou 

{Rj(u + ), Xj{u+) > s} U Mw+) - y(u_)} U {i2j-(u_), Xj(u-) < s} is a basis of M n . (2.21) 

The Majda stability conditions (2.20)-(2.21) are stronger than Lax entropy inequalities, and are satisfied 
by all Lax shocks in Systems (1.1) and (1.3) (see [38].) Majda's condition for contact discontinuities (of 
family k) is the following: 

{rj(u+), j < k} U {u + - u_} U {rj(u+), j > k} is a basis of E". (2.22) 

This condition is satisfied by any non trivial contact discontinuity in Systems (1.1) and (1.3). 

2.2 Interactions of weak waves, permutations of characteristic families and 
cancellation waves 

In this section, we consider the estimates for interaction of weak waves, that is, waves that are small. 
We begin by recalling the celebrated Glimm estimates. 

Proposition 2.1 ([24] & [18], Theorem 9.9.1). We assume that System (1.8) is strictly hyperbolic and 
satisfies (2.6). Consider (iti, 142,113) € £l 3 , and suppose 

u 2 = T(a',Ui), u 3 = T(cr",w 2 ), and u 3 =T(a, ui). 

Then 

J>, -ai- o'Dn = Y.aWiVi^} + 0(\o-'\\o-"\l\a'\ + \a"\}). (2.23) 

i j<i 

Moreover the "O" is uniform as u\, u 2 and U3 belong to a compact set o/Q. 

By the uniformity of the "0" as u±, u 2 and u 3 belong to a compact set of f2, we mean that, for some 
constant C > depending only on the compact K C where u\, u 2 and u 3 are chosen, one has 



5>*-^-<tf)ri-£^ , fa> r . 



3. 

j<i 



<C\a'\W"\[\a'\ + \a"\]. 



Note that the "— " sign in the right hand side in [18, Theorem 9.9.1] comes from its convention (7.2.15) 
on the Lie bracket: [r^Tj-] = (rj ■ V)r-j — (r^ • V)r^; here we prefer (as in [24] for instance) the convention 

[n,rj] = (n ■ V)rj - (rj • V)r,. 
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Remark 2.1. The point where r.i and [ri,rj\ are evaluated (among u\, 112 and u%) has no importance 
since the difference can be included in the term O(\o~'\\a"\ [|cr'| + |c"|]). 

Now an essential remark for the analysis developed here, is that it has no importance in Proposi- 
tion 2.1, that the characteristic families (Ai,rj) and the Lax curves T t are sorted in increasing order of 
the characteristic speed. This ordering of characteristic speeds only matters when translating the relation 
u u r = T n (a n , •) o ■ ■ • o Ti(o"i, ^ui" into an actual Riemann problem "find an self-similar entropy solution 
of the system with initial data u(0, x) = ui for x < and m(0, x) = u r for x > 0" . Incidentally, it is not 
important either that we use the usual wave curves T; rather than the wave curves from the right Tj 
(which is clear with our parameterization, see (2.15)) or the rarefaction curves IZi. 

A consequence of this is that one can permute the characteristic families, replace some Ti by Tt or 
IZi, and still get a result in terms of compositions of !&(•), Yfe(-) or IZk(-). Let us state this precisely. 
Let S n the n-th symmetric group. Given a permutation tt G S n and £ G {— 1, 0, 1}™, we let 

T^(«7, u) = T^> (a <n) , ■) o • • • o T$~™ (cr ff(1)l ■)«, o" € K n , u G fi, 
where we denote 

One can locally solve the "(7r, ^-swapped" Riemann problem exactly as in the classical case tt = id, 
£ = (1, . . . , 1) (this case corresponds also to tt : k 1— > n — fc, ^ = (— 1, . . . , —1)): given «_ and m + in O, 
close enough (at a distance which is uniform as U- lies in a compact of fi), there exists c G K™ such that 

m+ =T w ^(cr,w_). 

Indeed, all the maps involved are C 2,1 -regular as in the classical case. The argument relying on the 
implicit function theorem can be used without change. We denote 

cr = E 1r,f («_,«+), o-j = Y,J' i (u-,u + ). 

Now one can follow the proof of Proposition 2.1 in [18, Theorem 9.9.1] to check that the ordering 
according to the characteristic speeds does not intervene, and that the replacement of Tj with Tj merely 
implies to put a — sign before r^ in the estimates (due to the parameterization (2.16) of Tj). We obtain 
the following. 

Proposition 2.2. Let n G S n and £ G { — 1,0, 1}". We assume that System (1.8) is strictly hyperbolic 
and satisfies (2.6). Consider (1*1,1*25^3) G £1, and suppose 



u-2 



= T 7r '«(a', Ml ), U3=T*>t(<T",u 2 ), and u a =T**(a,u 1 ). 



Then 



J^in-rt-On = J2 (-l) J - 1 ^ w ^- 1 ^ 0) <(i)^a)[^(i),r w y ) ]+O(klk'|[K| + K'|]). (2.24) 

i 7r(i)<7r(i) 

Moreover the "O" is uniform as U\, U2 and U3 belong to a compact set of VI. 

When £ = (0, . . . , 0), this is the classical formula for permutations of flows of vector fields. 

An immediate corollary is that many waves conserve their nature (shock/rarefaction, increasing/decreasing 
contact discontinuity, compression wave/rarefaction wave) across an "interaction" . Actually, concerning 
Systems (1.1) and (1.3), one knows now in great details the result of the interaction of waves with large 
size: see in particular Chen and Hsiao [9] and Chen, Enders and Jenssen [10]. One of the issues in these 
papers is the possible appearance of vacuum, which is avoided here. 

Corollary 2.1. There is some n > uniform as u\ belongs to a compact set of VI, for which if 

u-2 = T[^(a' t ,ui), u 3 = Tp£(o-",un), and u 3 = T n, ^(a,Ui). 

with max(|cr^|, |cr'-'|) < n, and if i ^ j (resp. if i = j and a[, a'' have the same sign), then o~i has the 
same sign as o~[ (resp. as a' t , a" ). 
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We can deduce from Proposition 2.2 the existence of cancellation waves. By cancellation wave, we 
mean here a simple wave (u r ,u r ), associated to two simple waves (ui,u m ) and (u m ,u r ) and designed in 
order that, in the outgoing Riemann problem (u;, u r ), the wave associated to the characteristic family k 
vanishes. Here this simple wave takes the form of a rarefaction or compression wave. Precisely we obtain 
the following. 

Corollary 2.2. Suppose that n > 3. We assume that System (1.8) is strictly hyperbolic and satisfies 
(2.6). Consider (ui, 11,2,113) € Q, and suppose 

u 2 = Ti(a>i,Ui), u 3 = Tj(J)j,U2), 

with cti and /3j small. Let k G {l,...,n} \ {i,j}- There exists jk £ K such that, denoting a = 
£(ui,ftk(7k,u 3 )), one has 

o-k = 0, 
and additionally 

\*i - on] + \a 3 - /3j\ = OOoiH&l), ik = -Oi0jik ■ [n, rj ] + O (jaiWft] [K| + |/?,|]) . (2.25) 

Moreover the "O" is uniform as u\, ui and U3 belong to compact sets ofQ. 

Wc represent the result of Corollary 2.2 in Figure 1, where the waves are represented as single 
discontinuities. There can also be outgoing waves of "uninvolved" families, though we did not represent 
them. The case i = j is included in the result. 



no outgoing fc-wave 



j -wave 




j -wave 



'^fc(7fe>'"3) 



fc-wave 



Figure 1: A cancellation fc-wave 

Remark 2.2. As for Corollary 2.1, we can deduce from Proposition 2.2 some information on the nature 
of the additional wave (us,Tk(u3)): assuming that Ik • [ r iy r j] does not vanish on some connected compact 
set K, offl, there is some k > such that i/max(|a!j|, |/3j|) < K and u\ € K,, then 7^ has the same sign 
as —ctifijlk ■ [ri,rj\. In the same way, reducing k if necessary, o~i (resp. o-j) has the same sign as Ui 
(resp. /3j). 

We introduce also another type of cancellation wave in a sideways framework, which makes an "in- 
coming" wave vanish. 

Corollary 2.3. We assume n — 3 and that System (1.8) is strictly hyperbolic and satisfies (2.6). 
Consider (ui, 1*2,1*3) € Q, and suppose 

u 2 = Hi(-ai,u 1 ), u 3 = T 2 (f32,u 2 ). 

Then for some a e R 3 , one has 

u 3 = fci(-(Ti,-) o T 3 (er 3 , •) o T 2 (er 2 , •) Ui, (2.26) 

and additionally 

\a x - ai| + |a 2 - /3 2 | = 0(\ ai \\fo\), a 3 = ai f3 2 £ 3 ■ [n,r 2 } + O (|ai||&|[|ai| + |/3 2 |]) . (2.27) 

Moreover the "O" is uniform as u\, u 2 and u 3 belong to compact sets ofQ. 
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See Figure 12(c) below for a graphic representation — there ui, u m and u r replace u\, u-i and M3. 

Remark 2.3. As before, assuming that £3 ■ [ri,r2] does not vanish on some connected compact set K. of 
Q, there is some k > such that ?/max(|ai|, I/S2I) < K and u\ G K,, then o~\, o~q, and 03 have the same 
sign as ct\, /3% and —olx^^z ' [ r ij r 2], respectively. 

Proof of Corollaries 2.2 and 2.3. Corollary 2.2 is obtained by using the permutation 

1 ... k — 1 k ... n — 1 n 



" ' 1 ... fc-1 fc + 1 



and the vector £ = (1, . . . , 1, 0). Note in particular that one has u-i = T^^(cr' , ui) with o~' k = Ski&i an d 
u 3 = T v ^(a", u 2 ) with cr^' = 6 k jPj- Then 



")fc 



: = -S^' e («l,«3), 



satisfies the properties. 

Corollary 2.3 is obtained by using the permutation 

/ 1 2 3 

n= {2 3 1 

and the vector ^ = (1, 1, 0). □ 

2.3 Strong discontinuities, Riemann problem and interaction estimates 

Now we present some material allowing to work with strong discontinuities (shocks or contact discon- 
tinuities) in BV solutions. The material that we present here is mainly extracted from Schochet [38]; 
we recall it for better readability and to be able to be more specific on some particular aspect (see 
Remark 2.6 below). 

The first point is that the Riemann problem is locally solvable near a Majda-stable shock or a Majda- 
stablc contact discontinuity. 

Proposition 2.3. We assume that the system (1.8) is strictly hyperbolic and satisfies (2.6), as well 
as (2.8) for linearly degenerate fields. Consider (u_,u+) which is either a Majda-stable shock or a 
Majda-stable contact discontinuity: 

u + = T(a,u_), a := (0, . . . ,0,5^,0, . . .,0). 

There exist two neighborhoods w_ and uj + of U- and u+ respectively, a neighborhood S of a, such that 
for any u_ G cj_, any u + G uj + , the Riemann problem (v,-,U+) is uniquely solvable with a wave vector 
in S . Moreover there is a constant C > such that for all u l _, u 2 _ in uj_, all u\, m^ in u+, if 

u\ = T n (o-l-)o...oT 1 (o- i 1 ,-)uL, i = l,2, 

for wave vectors (<Tj)j=i.„ n an d {oj)j=i...n ^ n S } then 

n 

E w* -°'\< c (\ u - - u -\ + 1«+ - «+D- ( 2 - 28 ) 

Proof of Proposition 2.3. As in Lax's proof in the case where all waves have small amplitude, this is a 
consequence of the inverse mapping theorem. To establish the first claim, it suffices to check that 

dT 

-pr—{a, U-) is non-singular. (2.29) 

OCT 



It is elementary to check that 



q t 1 Ti(u + ) for i > k 

arfo«-)=i dT k (a k ,-) (2.30) 

CCT.J — rAu-) for 1 < k. 
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• 



Now let us begin with the case of a shock. To compute ^^-{a^uJ) and ^-(<r, u_) = ^r-(a,U-), one 
differentiates the Rankine-Hugoniot relation (2.9) to get 

dT k ds k (df _dtp\~ 1 r , 

d^^ U - ) = d^\d^- S d^j + ^' (2 ' 31) 

( u -) H-5T-* 1ST? 1^- s 7^r +b(«) J®-5T- • ( 2 ' 32 ) 



Here we used the notations 

[<p(u)] = <p{u+) - >p(u-), (2.33) 

and Sfc = Sfc(crfc,u_) = s(u—,Tk((Tk,u—)) and the index +/— means that the function has to be computed 
at u + /«_. Recall from (2.2) and (2.20) that the matrices < g£ — s ^ > are non-singular. We deduce 
that 

{(Xi(u + ) — s)Ri(u + ) for i > k, 
£[*>] fot . = , (,34, 

(A»(«_) - s)Ri(u-) + (n ■ V u s fe ) [<p(u)] for i < k. 

The assertion (2.29) is now a direct consequence of (2.21): by the inverse mapping theorem, the mapping 

A : (u_,a) !->•(«_, T(0-,u_)). (2.35) 

is locally invcrtible near (m_ , er) . 

• For what concerns the case of a contact discontinuity, one has 



q t j n(u + ) for i > k 

Tj(u_) for i < k. 



du 



Using (2.12), we infer 



{ri(u+) for i > k, 
r k for i = k, (2.36) 

ri(u-) for i < k, 

and one can conclude using (2.22). 

Now as before we denote £ the above mentioned inverse of A, that is, the mapping (u_,u+) G 
cj_ x uj + — > 5 defined by 

<T = £(u_,U + ) <^=> M+=T(<7, M_). 

Then estimate (2.28) is just a consequence of the Lipschitz character of S. □ 

Remark 2.4. Shrinking w_ and lo+ if necessary, we can assume that any simple wave with endpoints 
in w_ and lo + is automatically a Majda-stable discontinuity of family k and that each Riemann problem 
between two states in w_ or between two states in uj + is solvable. 

Remark 2.5. Equivalent formulas can be derived for wave curves from the right. 

The next point is to give interaction estimates for such large discontinuities. 

Lemma 2.1. We assume that the system (1.8) is strictly hyperbolic and satisfies (2.6), as well as (2.8) 
for linearly degenerate fields. Consider (u—,u+) G uj- x lj + which is either a Majda-stable shock or a 
Majda-stable contact discontinuity: 

«+=T(ff,«_), <f:=(0,...,0,cf fcl 0,...,0). (2.37) 
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1. (Interaction on the right) Let u r in lj + , and introduce a' :— £(it+,w r ) and a :— S(u_,u r ). Then 

° = °+(j^) ' rt°'Mu + )\+0{\u r -u + \ 2 ). (2.38) 

2. (Interaction on the left) Let ui in w—, and introduce a' : — £(u;,u_) and a := S(u;,u+). Then 

f n \ 

£>«»•<(«-) )+0{\u_-u l \ 2 ). (2.39) 



dT\ 1 /ar 



da I \du , . 

Proof of Lemma 2.1. This is a direct consequence of Taylor's formula for S, whose partial derivatives 
are easily computed from (2.35), and of Lax's estimates 



u r — u + 



= ^°-[n{u+)+0{\u r -u+| 2 ), or u- -ui = -^2ain(u-)+0(\u- -ui\ 2 ). 



□ 



Remark 2.6. We notice that is has no importance that the actual wave curves Tj are used in the 
Riemann problem a' := £(^2,143). We could replace them by TZi without change or by Tj with a mere 
change of sign on r,. In other words, we could consider a' := T^'^(u2,u^) instead and obtain the same 
result on a :— T^(ui,Ua), up to a change of sign if £j = — 1. The same is valid for an interaction on the 
left. 

Remark 2.7. Note that from (2.17), one infers 

'dTY 1 fdTY 1 f dT 



da J \da J \ du + J 

Hence the two formulas (2.38) and (2.39) are "inverted" when one replaces T with T. 

Now we distinguish the cases of a shock and of a contact discontinuity. 

Lemma 2.2. In the situation where the k-th family is GNL and that (w_,u + ) is a Majda-stable shock, 
we have a — a + a, with: 

• Case 1, Interaction on the right: 

{fu-' S %} + {^) a ^ a[{MU+) m(U+) + ° il<7 ' n {2Al) 

• Case 2, Interaction on the left: 



9/ __chp 
du du 



} (£) a = E^{(A,(«_) - t)Ri{u.) + (r<(u_) ■ Vs k )Mu)}} + 0(\af). (2.42) 



Moreover the "O" are uniform as ui,U- and u + ,u r belong to compact sets of w_ andu>+. 

Proof of Lemma 2.2. This is a direct consequence of Lemma 2.1 and of (2.31)-(2.32). □ 

Lemma 2.3. In the situation where the k-th family is LD, that (2.8) is satisfied and that (u-,u+) is a 
Majda-stable contact discontinuity, we have a = a + a, with in the case of an interaction on the right, 
respectively on the left: 



a = Y J <r l {u + )+0{\a'f), resp. — a = £<#■<(«_) + 0(\af). 



Moreover the "O" are uniform as u;,u_ and u + ,u r belong to compact sets of w_ andui+. 

Proof of Lemma 2.3. This is a consequence of Lemma 2.1 and of (2.36). □ 

Note that the matrices appearing in Lemmas 2.2 and 2.3 have been computed in the proof of Propo- 
sition 2.3, and that one can use T rather than T, which in some situations can simplify the writing. 

15 



3 About systems (1.1) and (1.3) 

3.1 Some characteristic elements of systems (1.1) and (1.3) 

3.1.1 Eulerian system (1.1) 

For the Eulerian system (1.1) one fixes: 

u= \v , «€!l = M + xIxK+ (3.1) 

ip(u) = pv I , /(«) = pv 2 + P , (3.2) 

V V^ 2 + P ) \ V/ 9 ^ + l Pv J 

so that (p maps primitive coordinates to conservative ones. We have 

1 \ / 1 \ 

dip(u) = | v p , dip-^u) = -v/p l/p , (3.3) 

2=V {j-l)pv 1/ V 2=V -(7-l)v 1/ 

u p 0\ / v p \ 

d/(u) = ( v 2 2pv 1 , (dvj)- x 4f(«) = v l/p . (3.4) 

V 1 ^ 3 |(7~ 1 )/°« 2 + 7 p 7 U / V ° l p v J 

The characteristic speeds (i.e. the eigenvalues of (dip)~ 1 df) are given by 

Ai = u - c, A 2 = v, A 3 = v + c, (3.5) 

with the speed of sound c given by 




The eigenvectors of {dip) l df(u) are given by: 



Ra = -, ^ ; I '•-'' 



(7 + l)c , ^ 

v ' ' \ (7 — ljuc — c 



2 - V- 2 



(3.6) 



9 / -P/ C \ ( 1 \ o / P/ c 

and those in terms of the ip variable (that is, the eigenvectors of df(dp)~ 1 , i.e. Ri = dip ■ r,): 

(7+1)c V h-l)vc + c 2 + V« 2 / V ^« 2 

1 



These eigenvectors r^ i = 1,3, satisfy in particular (2.7). The corresponding left eigenvectors of 
(dp)~ 1 df(u) are 

*- ^ (<4 -£) • *-('■ »• 4) ■*- ^ (°- 1- ^) ■ <-> 

It will be useful to extend the definition of the shock speed to rarefactions/compression waves. If 
u 2 =Tli((Ti,ui), Ui = (pt,Vi,Pi), i = 1,2, we set 

P2«2-P1W1 J Q CTl A l (^ i (cr,Ui))i? l 1 (^(c r ,Wl))rfo- 

s(mi,-u 2 ):= = rCTa Dl ,„ -. ^~, , 3.9) 

P2-P1 J R\{Hi(<j,u\))da 
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where R\ stands for the first coordinate of Ri . This obviously gives the shock speed for actual shocks as 
well. It is clear that this shock speed also satisfies (2.10), with a O uniform on compacts of fi and that 

s(ui,u 2 ) e [\i(ui),\i(u 2 )} U [Xi(u2),Xi(ui)]. (3.10) 

Note finally that in the coordinates given by (3.1), the 2-contact discontinuities are given by 

U + = M_ + OT2, CTgl, 

and in particular v and P are preserved across the discontinuity. 

3.1.2 Lagrangian system (1.3) 

For the Lagrangian system (1.3) one fixes: 

v), (p(u)=\ v \ f(u)=\ P |, «e!l = R + xlxI + . (3.11) 

PJ V^ + W \Pv 



We have 



1 \ / 1 

dip(u) = | 1 , dip- l {u) = 1 ) , (3.12) 

A v =Fl/ V-f "(7-1)? 

-1 \ / -1 

df{u) = | 1 , {d V )- l df{u) = 1 ) . (3.13) 

P v \o^O 



The characteristic speeds are given by 



hyP c yP c 

\i=-\ /— = — , A 2 = 0, A 3 = a/ J - = -. (3.14) 



T T \ T T 

with c given again by (3.6), that is 

c= a/7-Pt 

The eigenvectors of {dip)~ l df (u) are given by: 



ri = ^| r J,„-l [U. ^1 / |. P.5, 



These eigenvectors r,-, i = 1,3 satisfy (2.7). The eigenvectors in terms of the <p variable (that is, the 
eigenvectors of df(d(p)~ 1 , i.e. Ri = dip ■ n) by: 

*-?M-,!J-*-(i)'*-o^U-)- (3M) 

The left eigenvectors of (dip)~ 1 df(u) are 

and the ones of df(d<f)~ 1 (i.e. Li — ti- (dip)^ 1 ): 

*-( 3 r i .-fr-i>£. 3 £). ( " 9) 

^ = ^-Al-(7-D^). (3.20) 

4r \ 7r c c J 
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Here we extend the shock speed to rarefactions/compression waves as follows: if u 2 — 7?.j(<7j,ui), u% — 
(Ti,Vi,Pi), i = 1,2, we set 

/ s v 2 -v 1 J^ \i(1li(a,ui))R}(lli(o;ui))dcT 

s(ui,u 2 ) := — = ai — , (3.21) 

t 2 ~t x J Q R\ (Jli {a, ui ) ) do- 

where R\ stands for the first coordinate of Ri. Of course, (2.10) and (3.10) apply here as well. 

Finally it will be useful to have the the shock/rarefaction curves described in the Lagrangian case 
(though in fact they coincide with the ones in the Eulerian case after change of variables.) They can be 
parameterized through the coefficient 

P+ 

x= pZ 

as follows (see e.g. [40, Section 18. B], [11, Section 6.4]): the shock curves are given by u+ = Si(x,uJ) 
with 

' P+=.xP_, 

T + /3 + X 

tZ = Jx + i ' (3.22) 

/ 2 1-x 



< 



_ 7(7-l)vWTT : 
with x > 1 and the + sign for the 1-shocks, with x < 1 and the — sign for the 3-shocks; we have put 

0:= 1+1. (3.23) 

7-1 

Note that 

[P] > 0, [t] < 0, [v] < across a 1-shock, [P] < 0, [t] > 0, [v] < across a 3-shock. (3.24) 
The corresponding shock speed is given by 



c- l + px c+ / P + x 

T ^VTTJ = T ^V(TT^- (3 ' 25) 



The rarefaction curves parameterized by x are given by u + = lZi(x:ii-) 



P+=xP_, 

— — lAy— 

T+ = X ' 'T- 



2c- 1 t 

V + = V- ± (x^ — 1). 

7-1 



(3.26) 



with x < 1 and the — sign for the 1-rarefactions, with x > 1 and the + sign for the 3-rarefactions; we 
have put 

(:=^1. (3.27) 
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Rarefactions of the non-isentropic Euler actually coincide with curves of the isentropic model; the physical 
entropy S is conserved along those curves. 

The family of 2-contacts discontinuities is simply described by u + = u_ + crr 2 , so that one has 

[P] = 0, [v] = across a 2-contact discontinuity. (3.28) 

3.1.3 Commutation of rarefaction and compression waves of families 1 and 3 

An important relation satisfied by both systems is that 

4-[n,r 3 ]=0. (3.29) 

This does not mean that r± and r 3 commute, but they satisfy the integrability condition of Frobenius' 
theorem. It follows that the curves IZi and IZ3 locally define a submanifold of K 3 of dimension 2, on 
which (0-1,(73) H> lZi(<7i, ■) o 7^3(0-3, -)u and (ctl, 03) H> 7^3(0-1, •) o TZ\(ai, -)u give local diffeomorphisms. 
A consequence is the following. 
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Lemma 3.1. Consider both Systems (1.1) and (1.3). Let ui £ £1 and o\ G K and (73 G K small. 
Let 7r e 53 anrf £ = (0,0,0). If u r = lZi(ai,-) o 7^3(173, -)ui or u r = TZ^a^,) o 7£i(<7i, -)ui, then 
Y12 ( u h u r) — 0, where £2 designates the second component ofT,™'^. 

In other words, the interaction of rarefaction/compression waves of families 1 and 3 does not generate 
a 2-contact discontinuity (as long as one considers the Riemann problem in terms of compression waves 
rather than in terms of shocks.) Note that for large rarefactions, one has to be careful about the possible 
appearance of the vacuum (see [10]). 

Another way to look at Lemma 3.1 is to notice that the physical entropy S is constant along the 
curves IZi, i = 1,3. This can be seen by a direct computation or relying on (1.6) and (1.7) which are 
deduced from (1.1) and (1.3) for regular solutions. Hence the submanifold mentioned above is a level 
surface of S (on which the solutions of the isentropic equations live, by the way). But it is obvious 
that following the wave curve Ti = 7?. 2 of the second family increases/decreases S, so there cannot be a 
non-trivial second component in E^'S 

3.1.4 Notations for the elementary waves 

We will use the same notations as [11] and [10] to describe elementary waves. The notations are as follows: 
S will designate a shock, R a rarefaction wave (either of the first or of the third characteristic family), 
J will designate a contact discontinuity (of the second family). We add C to designate a compression 

wave. Waves of the first family will be more precisely described as S , R and C, those of the third family 

S , R and C. We will distinguish between the contact discontinuities satisfying t_ < r + where t_ (resp. 

< 
r+) is the specific volume on the left (resp. on the right) of the discontinuity, which we denote by J, and 

> 
those for which t_ > t + , denoted J. We underline that we use this notation including for the system in 

< 
Eulerian coordinates for which we rather use p as an unknown; in particular a J satisfies p_ > p + and 

it corresponds to u + — T2{o~2, uJ) with a 2 < 0. 

In the figures, in order to emphasize the waves that we consider strong, we will put in this case these 

letters in blackboard bold style. In those figures, we may also use the letter A to designate "artificial" 

waves (see below). 

3.2 Some coefficients for interactions with strong discontinuities 

Here, we compute several coefficients allowing to estimate the strength (and the nature) of outgoing waves 
for some particular strong discontinuity/weak waves interactions, using the tools of Subsection 2.3. In 
the case of the Eulerian system, we are particularly interested in the interaction of a small 3-wave with 
a strong 2-contact discontinuity. In the case of Lagrangian coordinates we are particularly interested in 
the interaction of a small 3-wave with a strong 3-shock (or the interaction of a small 1-wave with a strong 
Tshock which can be deduced from the latter through the change of variable x < — > —x associated with 
(r.v.P)*— ►(t,-«,P)). 

Notation. The coefficients that we will introduce connect the strength er' of a weak wave of family j, 
interacting with a strong wave of family k (of strength a^), with the strength oi of the outgoing wave of 
family i. This will be written as 



' |2\ 



<7i = 5 ik a k + {a l k] or a)^ o' d + 0(\a 

where akj (respectively aj_k) means that the j-wave interacts with the strong fc-wave from the right 
(resp. left). The coefficient 8ik is Kronecker's symbol. For instance in (3.31) below, the coefficient 
a\ 2 appears when computing the strength of the outgoing wave of the i-th family as one considers the 
interaction of a weak wave of the third family with a strong wave of the second family, the weak wave 
being the left one. 

3.2.1 Eulerian case 

Here we consider, in the case of System (1.1), the interaction of a 2-contact discontinuity (ii-,u + ) 
(considered as strong) with a wave of the third family, situated on its left. Of course the case of a strong 
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2-contact discontinuity interacting with a 1-wave on its right is similar (and obtained via the change of 
variable x < — > —x.) We prove the following. 

Proposition 3.1. In System (1.1), let u_ = (j>_,V—, P-) and u + = (p, ,u+,P+) two states in Q = 
R + x R x R + joined through a 2 -contact- discontinuity: 

u+ = T(a,U-), ct = (0, ct 2 ,0), a 2 ^ 0. 

We consider w_ and w + as m Proposition 2.3. Let ui,U— in w_ and n + m cj + satisfying 

u+=T(a,u-), &= (0, a-2,0), <r 2 7^ 0, 
w; = T(cr',u_), cr' = (0, 0, 0-3). 



TTien one /las 






with 








o\ 


= ^20-3 + 0(0 


where 








n- 1 - 


_ y/P+~ ^/P 1 



S(n;,M H 



cr + cr, 



o-2 = al g o-' 3 + 0(a' 3 2 ), 



ip- 



2 v^ 



*3,2 



*3,2 



*3,2 



(3.30) 
(3.31) 

(3.32) 



Moreover the "O" are uniform as u;,u„ and u + belong to compact sets of u)_ andw+. 
Note in particular that a 3 2 > (resp. aj 2 < 0) when ,o_ < p + (resp. when p_ > p + ). 

Proof of Proposition 3.1. First, it is straightforward to see that Majda's stability condition (2.22) is 
satisfied by (u_, u + ). Now according to Lemma 2.3 and to (2.36), one has S(u;, u + ) = a + cr, with 

ain(u-) + cr 2 r 2 + a 3 r 3 (u + ) = a 3 r 3 (u-) + C(|cr 3 | 2 ). 

We consider the matrix 

^i(«_)-n(«_) £i(u-)-r 2 &(«_)• r a (ti+)\ (l 4(w_) • r 3 (n+) N 

<a(«-)-ri(tt_) 4(u_)-r 2 * 2 (u_) • r a (u+) = 1 *a(u_) • r 3 (u+) 

^a(«-)-ri(u_) 4(u-)-r 2 * 8 (u_) ■ r 3 (u+)J \0 * 3 (u_) • r 3 (u+), 

so that one has Per = a' + 0(|cr 3 | 2 ). Inverting P, we finally end up with (3.31) with 
1 h(u-) ■ r 3 (u+) 2 £ 2 (u-) ■ r 3 (u+) 3 1 

3,2 h(u_)-r 3 (u + y a3 ' 2 - 4(tt-)T 3 (u+)' 3 ' ? 4(tt-)T 3 (u + )' 

the denominator being always positive. Computing these coefficients leads to 



P :-- 



a 



3,2 



p+c+ - p_c_ 
P-C- + P+c+' 



l 3,2 



2p_p + (c 2 h -c 2 _) 

c+(p_c_ + p+c+y 



a 



2p_c_ 



3,2 



(p_c_ +/)+c+)' 



But since the pressure is constant across a 2-contact discontinuity, these formulae simplify to (3.32). □ 

Remark 3.1. The situation where a 1-wave interacts with a 2-discontinuity from the right is exactly 
symmetric. Note that in the symmetry x < — > —x, a wave 1-wave (resp. 2-wave, 3-wave) is transformed 
into a i-wave (resp. 2-wave, 1-wave) with the same (resp. opposite, same) strength. Hence one gets the 
same result as Proposition 3.1 with 

2^/P+ 2 2JJZ{JJZ - JJ^) 



"2,1 



"2,1 



"2,1 



Ip- 



(3.33) 



VP- + VP+ " VP- + VP+ " VP- + VP+ 

so a 2 1 > (resp. a 2 ± < 0) when p_ > p + (resp. when p_ < p+). 

Remark 3.2. Shrinking U- and uj + if necessary, we can ensure that, in the case of the interaction of 
a 3-wave on the left (resp. a 1-wave on the right) of the strong 2-discontinuity, o~\ has the same sign as 
a 3 2 cr 3 (resp. a\ x a[) and o~ 3 has the same sign as a 3 2 a' 3 (resp. a 2 1 a[). 

Remark 3.3. Using the notation f — 0(<?) to express that there exist c, C > (possibly depending on 
U-) such that, for small values of the variables, 

cg<f< Cg, (3.34) 

we have, for p_ > p + that a\ 2 = 0(— <t 2 ), oc\ 2 = Q(a 2 ) and a 3 2 = O(f). 
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3.2.2 Lagrangian case 

Here we consider, in the case of System (1.3), the interaction of a 1-shock (u_, u + ) (considered as strong) 
with a wave of the first family, situated on its right. Again, we are interested in estimating the resulting 
outgoing waves using the tools of Section 2.3. Of course the case of a strong 3-shock interacting with a 
3- wave on its left is again obtained via the change of variable x < — > —x, (t, v, P) < — > (t, — v, P). We 
prove the following. 

Proposition 3.2. In System (1.3), let 11- = (t_,TJ_,P_) and u + = (T + ,v+, P + ) two states in tt — 
R + x R x R + joined through a 1-shock: 

u + =T(a,u-), ct = (cti,0, 0), a\ < 0. 

We consider w_ and w+ as in Proposition 2.3. Let u_ in w_ and u + ,u r in lu + satisfying 

u + =T(<t,m_), (7=(d-i,0,0), di < 0, (3.35) 

u r --T(a\u + ), a 1 = K, 0, 0). 

Then one has 

£(u_, u r ) = a + a, (3.36) 

with 

a 1 = al l a' 1 +0(a[ 2 ), <r a = 0^ 1 a{ + O(<j^), a 3 = a\ A a[ + 0(a[ 2 ), (3.37) 

where, denoting s — s(u-,u + ), 

1 1 2 _ L 2 {u+) ■ [<p(u)] Xi(u+)-s 3 _ L 3 (u+) • [<p(u)] \i(u+) - 8 



^L 1 {u+)-[<p(u)]' l ' x Li(u+) ■ [<p(u)] s ' ^ L 1 (u+) ■ [<p(u)] X 3 (u+) - 



•s 



(3.38) 
Moreover the "O" are uniform as u_ and u+, u r belong to compact sets of W- and w+. 

Proof of Proposition 3.2. It is straightforward to see that Majda's stability condition (2.20)-(2.21) is 
satisfied by any 1-shock (u_,u_|_). We use Lemma 2.2 and (2.34); we have 

d j^-fo>(u)] - sa 2 R 2 {u+) + a 3 {X 3 {u+) - s)R 3 {u+) = a[(X 1 (u+) - s)R 1 (u+) + 0(\a[\ 2 ). 

We consider the matrix (using again the notation (2.33)) 

/^£i(u+)-[v(»)] -«ii(u+) • R 2 (u+) (A 3 (u + )-«)Li(« + )-fl8(«+^ 
P = & M«+) ' [¥>(«)] ~sL 2 (u + ) ■ R 2 (u + ) (X 3 (u + ) - s)L 2 (u+) ■ R 3 (u+) 
\m L z( u +) ■ b( u )] -sL 3 (w+)-^2(w+) (As(u+) - a)I 3 (tt+) • Ra(u+), 

^£i(«+)-b(«)] 

^£2(1*+) • [?(«)] -s 

^i 3 («+) ■ b(«)] A 3 ( U+ ) - 

Inverting P, we obtain (3.37) with the coefficients given in (3.38). □ 

Now we focus on the coefficients a\ 1 and a\ 1 . One has the following result. 

Lemma 3.2. Any 1-shock (3.35) satisfies 

al A > 0. (3.39) 

Moreover, for 7 < |, any 1-shock (3.35) satisfies 

al A < 0. (3.40) 

Shrinking u;_ and lj + if necessary, these coefficients are uniformly strictly separated from zero, and in 
(3.37), Ok has the same sign as a\ x a[, k = 2,3. 
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Proof of Lemma 3. 2. First one has clearly 

Ai(u + ) — s < 0, s < and A 3 (u + ) — s > 0, 

since the Lax inequalities are valid on the whole shock curve (see (3.25)). Now let us determine the signs 
of the various Li(u+)-[ip(u)]; actually it will be a bit simpler to work with sLi(u+)-[tp(u)] = Li{u + )-[f{u)\. 
Using the Rankine-Hugoniot relation, we find that 

[/(«)] = [ [P]\ = M ( » 

\[Pv]J \P- + sv +/ 

Now one computes sLi(u + ) ■ [(p(u)] as follows: 

Note that each term inside the parentheses is negative. 
Concerning sL 2 (m+) • [¥>(")], one has: 

sL 2 (u+) ■ Mu)} = -l— L ( v+ - v _)-( 1 -l)^ r (P + -P_) + l^ ) 1 (P +v+ -P_ v _) 

7 lP+ lP+ 

= 1^ -„)(£-! 

7 V p + 

Using (3.24), we deduce (3.39). 

The factor L^(u + ) ■ \<p(u)\ is the one sensitive to 7. One has 

sL 3 (u + ) ■ [<p(u)] = 1 -±l[ v ](^ + s + ^lp_ 
4r + V7^+ c + 

(7 + 1)04..,/. ,,P- 



Now we use the representation (3.22) of the 1-shock curve. We find that 



(7+l)c+ r 1 (7-1 1 I 13 + x 

rj 



sL 3 (u + ) ■ [ 9 {u)} = -^- M _ + 1 - 7 y ( IT ^ 



To determine the sign of the last factor, we parameterize the function by t = l/x and hence consider 
only t G (0, 1). The function 

ft(t):=( 7 -l)t+l-V '" " 



, (1 + /3) 
vanishes at 1; its derivative 

ft'(£) := (7-I) ^— , * - , 

is increasing and 

h'(l)=-y-l + l±±<0*=* 7 <| 

It follows that for 7 < ~, ft,' < on (0, 1) and consequently ft > on (0, 1). Hence we find that (3.40) 
holds for 7 < |. 

Finally, we obtain a negative upper bound for the coefficients a\ 1 and a\ 1 by choosing small compact 
neighborhoods w_ and w+ of w_ and w + inside of w_ and oj + such that any 1-shock from w_ to w + 
satisfies 

Ai(u+) — S < — k < 0, s < — k < and A 3 (u + ) — 3 > « > 0. 

n 
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Remark 3.4. This implies in particular that the interaction of a strong 1-shock with a small 1-shock 
generates a rarefaction wave in the third characteristic family, a fact that is well known, see e.g. [40, 



5 
3> 

complex, see Chen, Endres and Jenssen [10] . 

Remark 3.5. With the notation of Remark 3.3, using (2.10), we can see following the lines above that 

5 " 



Theorem 18.8]. For 7 > |, this generates a shock, but the interaction of two strong shocks can be more 



a \ x = 9(1), a? ! = ©(of) and (for 7 < §J a\ x = ©(-of). 



We will be also interested in the result of the interaction of such a strong 1-shock with a weak simple 
wave (of family 1, 2 or 3) on its left. 

Proposition 3.3. In System (1.3), let u_ and u + two states in Q = R + x M. x R+ joined through a 
1-shock: 

u + =T(a,u^), (j = (<7i, 0, 0), <7i < 0. 

We consider uj^ andto + as in Proposition 2.3. Let i G {1,2,3}. Letu^,ui inuj^ and u+ inuj + satisfying 

u+ = T(ct,u_), <t = (0-1, 0, 0), <7i<0, (3.41) 

m:=T 4 (^,«_). (3.42) 



Then one has 
with 



£(u;,u-i_) = o + cr, (3.43) 

C7 1 = al, 1 ^ + 0(c7f) ) os = a? 1 oi + 0( i a ) ) <t 3 = 0?^ + 0(of ), (3.44) 



w/ie 



af 5l > 0, a? x < («/7<-), a| 4 > 0, a| x < 0, a^ < awd af^ > 0. (3.45) 



TTie coefficients af 2 are bounded and bounded away from zero uniformly and the "O" are uniform as 
Ui , w_ and u + belong to compact sets of w_ and uj + . Moreover, shrinking lo^ and lo + if necessary, erj, 
has the same sign as oq x a^ in (3.44). 

Remark 3.6. We can deduce as before the equivalent for weak waves interacting on the right of a strong 
3-shock through the change of variables x < — > —x. Recall that this transforms a 1-wave (resp. 2-wave, 
3-wave) into a 3-wave (resp. 2-wave, 1-wave) with the same (resp. opposite, same) strength. 

Proof of Proposition 3.3. We follow the same lines as before, applying Lemma 2.2 and the formula (2.34). 
We get that 

B 3 

*!■£ b(u)]+5>j(Aj(u+)-«)i^(u+) = -^(A,( U _)-s){i? i (u_)+(r i ( U _)-V„ S ) [tp(u)] }+0(K| 2 ), 

1 J =2 

The matrix whose columns is formed by [y(u)], R 2 (u + ) and R 3 (u + ) is invertible. Hence we find that 
(3.44) holds, and in particular we can compute the coefficients 

2 _ K(u-) - s det(Rj(u_),R 3 (u + ), [ip]) 3 _ Aj(u_) - s dct (R l (u_),R 2 (u + ),[p]) 

a ^-~ s det{R2{u+),R 3 (u+),[<p}) " m< ■ a *±- \ 3 (u + )-sdct(R 2 (u + ),R 3 (u + ),[<p})- 

The quotients »(^-)~ s anc [ x (u~ V^-l are c l ear ly positive for all i; let us determine the signs of the 
determinants. Let us remark that using the Rankine-Hugoniot relation, we can replace [p] with [/] in 
these determinants. To simplify the writing, we will compute the determinants with the following vectors 
Ri and -R3 instead of R\ and R 3 : 

- _ ( 7 +i)c / ; \ . . ( 7 +i)c 
R1 -~to*~ R1 -[_p\e. )' R3 -^^ R3 ~ 

Since we are only interested in the signs of the determinants and since the factor ^ 2 is positive, 
bounded and bounded away from on compacts subsets of f2, this replacement is harmless. 
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We obtain: 

dct(R 2 (u + ),R 3 (u + ), [/]) = C -±[Pv] (p + + ?±?±) [P] -^-[P] + -^±-^[v] 

= ^(P- + ^)[v]-^[P]<0, 

r + \ 7 — 1 / 7 — 1 

since both [v] < and [P] > across a 1-shock. 
• Weak wave of family i = 1. We compute: 

del (/?,(„ )./?,«/+).[./■]) = -['-^— - : U ](^^-^±)-[P] 2 -[P](^±-^^)+[P,]^'- ' 



c_c + 


M^ P+C " 


P_ c+ 


c_c+ 


_ rpi 2 ^ n 


^+ 



Hence we deduce that a\ 1 > 0. 

det(R 1 (u-),R2(u + ),[f]) = -—([Pv] + ^[v])-[P](^ + P--^ 

r_ V 7— 1 / V7— 1 t_ 



r_ 7 — 1 V7— 1 

We use [P] = s[v] and the formulae (3.22) and (3.25) to obtain, with x = P+/P- > 1: 



det(R l (u-),R 2 (u + ),[f}) = -^^! : /' i+; -'-/ l ' 



7— 1 y 1 + /3 \7 — 1 x 
We define 

and observe that h(l) = and that 

1 ( P 2 , P. 



h ' {x) VT+PVT+pExi { 2( 7 -l) 1 + 2 X + 1 

In particular ft'(l) = "rig ^li^ ^ s ne g a tive whenever 7 < |, and one checks that in that case 
h! is negative on (l,+oo). Recalling that [v] < 0, we deduce that, provided that 7 < |, one has 
dct (Ri(u_) , R 2 (u + ) , [/]) < 0, and hence that a\ 1 < 0. 

• Weak wave of family i = 2. Next: 

det(R 2 (u-),R 3 (u + ),[f}) = C ^[Pv}-(p + + C -^±)[P} + ^—(^±[v}-[P] 

= ^-^[v]-(p + + ^-)[P]<0, 
7 — 1T+ V 7 — 1 / 

where we used (3.24). Hence a 2 1 > 0- Now: 



[PI 2 

det(iZ 2 (ti_) J i$ 2 (u+) ) W) = — LJ t<0. 

7-1 



It follows that a| 1 < 0. 
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• Weak wave of family i = 3. We have: 

det(R 3 (u_),R 3 (u + Uf]) = -M 2 ^-[<^-^W] 2 + [P](^-— )+[Pv](±* 

T-T+ V T - T + 

Using (3.25), we see that, with x > 1, 

It follows that a\i < 0. Finally: 

det( J R 3 ( W -),i? 2 ( U+ ), [/]) = -"— ([Pv] + ^jM) + [^ (^1 +P - 

reasoning as before. Hence a| 1 > (recall that s < 0). 

This ends the proof of Proposition 3.3. □ 

Remark 3.7. .Here using the notation of Remark 3.3 and following the lines above, we deduce that 
"1,1 = ©(^i), "1,1 = ©(-^ 2 ) 0-7 < §;, "1,1 = 9(1), a|,i - e(-^i), a|,i - 0(-*i) and a 3 3A - 6(1). 

3.3 Additional cancellation waves and correction waves 

In this section, we introduce additional cancellation waves and what we will call correction waves, relying 
on the strong discontinuities for Systems (1.1) and (1.3) described above. These waves will be compression 
waves. 

3.3.1 Eulerian case 

As in Paragraph 3.2.1, the strong wave that we consider in the case of System (1.1) is a 2-contact 

discontinuity. We look for compression waves of the third family, on the left of the 2-contact discontinuity, 

that cancel the effect of a 1-shock interacting on its right, in the sense that there is not outgoing 1-wave 

after the interaction; see Figure 2 (where the compression wave is represented as a fan of discontinuities 

< 
focusing to a point). The 2-contact discontinuity that we consider is a J wave (that is, for which r_ < 7+) 

> 
rather than a J one. 



C 




Figure 2: A 3-comprcssion wave acting as a cancellation wave 
Precisely, we establish the following proposition. 
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Proposition 3.4. In System (1.1), consider a 2-contact discontinuity (u—,U+) as in Proposition 3.1 
with a 2 < 0. We consider w_ and uj + as in Proposition 2.3. Let m in U— and u m , u r in u + satisfying 

u m =T(a,ui), cr= (0, <t 2 ,0), a 2 < 0, 
u r = T(0,u m ), = (0l, 0,0), /?! < 0, 

with |/3i| small. There exists 73 < such that, if ui :— 72.3 (— 73, u i), denoting a = T,(ui,u r ), one has 

o-i=0, 

and additionally 

73 =-^01 + 0(101 1 2 ), t7 3 =(al il -42^r i ) /9i+0(|ft| 2 ). (3.46) 

a 3,2 \ a 3,2/ 

Moreover the admissible size of \/3i\ and the "O" are uniform as ui, u m and u r belong to compact sets 
of CO- and lu + . 

By the uniformity of the admissible size of \Pi\, we mean that ui and u m being in fixed compacts of 
o;_ and uj + , there is a j3 > such that the property is valid whenever |/?i| < /?. 

Proof of Proposition 3.4- Given m and <7 2 , we consider the mapping : 

G : (/3i, 73 ) G (~£,e) 2 -> (ft, £i(ft 3 (-73, «l).Ti(|9i,r 2 (a a ,«j)))) , 

where as before Si denotes the first component of £ and e is a small positive number. It is clear that G 
is C 2 and its differential at (0, 0) is given by 

dG(o,o)=n ° 

\ a 2,l a 3,2 

From (3.32) we know that a\ 2 ^ and in fact, since a 2 < (so that p/ > p m ), we deduce that a\ 2 < 0. 
Even, one can have a negative upper bound for this coefficient, shrinking w_ and w+ if necessary. It 
follows that dG(0, 0) is invertible, of inverse 

dG(0,0)- 1 =(_ a | !i i j. (3.47) 

Hence the existence of 73 for ft small is the consequence of the inverse mapping theorem (and one can 
bound from below the size of ft for which this is possible in terms of |jdG(0,0) _1 || and ||G||cra). The 
first estimate in (3.46) follows from (3.47). The second estimate on 

(73 = E 3 (ft3(-73,tti)>Ti(/3i,T 2 (a- 2 ,u,))), (3.48) 

is then a first order Taylor expansion. That 73 < comes from the computation of the coefficient 



*2,1 



2 v/^> 



a 3,2 VPm - y/Pl 

which is positive since, is the case considered here, one has pi > p m . □ 

Remark 3.8. The sign of a 3 for |ft| small is given by (3.46); the coefficient can be computed as 



a 



3 "2,1 



1 I nr~ r^\1 



{y/Pm ~ yfplY + 4y/piy/p 



"2,1 _ "3,2 ~T 

"3,2 Pl-Pm 

which is positive so 03 < 0. However this is not essential in the construction and we will not use this 
fact. 
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3.3.2 Lagrangian case 

For the system in Lagrangian coordinates, we consider not only cancellation waves, but also compression 
waves which do not cancel one of the outgoing waves, but rather force the outgoing waves to have a 
prescribed sign. Hence we refer in that case to these compression waves as correction waves. 
We prove the following result (see Figure 3). 

Proposition 3.5. We consider System (1.3). Let us be given a 1-shock (u_,u + ), u + = T\(W\,u_) with 
<7i < 0. We consider w_ and lu + as in Proposition 2.3. Let u m in u;_ and u r in uj + satisfying 

u r = T(a,u m ), <r = (cri, 0, 0), a x < 0, 

and consider ui in uj_ such that ui = Y(/3, u m ) with j3 a simple wave of the form (/3i , 0, 0), (0, 02, 0) or 
(0, 0, /3s), with \/3i\ small. 

Then there exists 7i < such that, if u r := TZi(^i,u r ), denoting a = T,(ui,u r ), one has 



a 2 < 0, ct 3 > 0, 



(3.49) 



| 7 i|=0(|ft|). (3.50) 

Moreover the admissible size of \(3i\ and the "O" in (3.50) can be taken uniform as ui, u m and u r belong 
to compact sets of uj- andu>+. 




(a) Correcting the effect of a 1- 
rarefaction 



(b) Inverting a J wave 



(c) Correcting the effect of a 3-shock 



Figure 3: Compression waves acting as correction waves 



Remark 3.9. The cases /3 = (/3 1; 0, 0), ft > 0, /? = (0, /3 2 , 0), (3 2 > and (3 = (0, 0, /3 3 ), ^ 3 < 0, are 
the non trivial cases (for the other ones, 71 = works). In these situations, we could make sure that 
one of (T2 or 03 actually vanishes. But it would not be systematically the same one, so this information 
is of no use to us in the construction. 



Proof of Proposition 3. 5. Let 



and let us prove that 



/./,, 



"1,1 


) 


«!,! 


"1,1 


a l,l 



7i := 



-ilH + l)|ft| 
works. We compute er^, k — 2 or 3, by 

o-fc = Gkifiul) '■= Sfc(Ti(/3 l ,u m ),^i(7i,Ti((Ti,u m )). 

As in the proof of Proposition 3.4, Gk is of class C 2 and its differential at (0, 0) is 

dG k (0, 0) = a k iS df3i + a^ cfy. 

Hence 

G fc (ft,7i) = a?iA + a^7i + 0{\^\ 2 + | 7 i| 2 ). 



(3.51) 
(3.52) 



(3.53) 
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Since a\ 1 is positive and a\ 1 is negative (see Lemma 3.2), one deduces that 

G 2 (A,7s) < -IAI +0{pj) and G 3 (A,7 3 ) > |AI - O{0f), 
which allows to conclude. 



(3.54) 

D 



We will also consider actual cancellation waves for System (1.3), as stated in the next proposition. 
We refer to Figure 4. The strong shock used here is chosen not too large. 

Proposition 3.6. In System (1.3), there exists k > such that the following holds. Consider a 3-shock 
(w_,m + ) given by u + — T 3 (a 3l u-) with a 3 £ (— k,0) . Let uj- and lo + as in Proposition 2.3. Let ui in 
<jj- and u m , u r in lo + satisfying 

u m = T(a,u{), a = (0, 0, 03), 0-3 < 0, 
u r =T(f3,u m ), 0= (0,02,0), 2 <O. 

Suppose, shrinking w_ and lo + if necessary, that 173 £ (— ^f ,0). Then for small \f3 2 \, there exists 73 < 
such that, if ui :— 72-3 (— 73, Ui), denoting a = S(w;,u r ), one Acts 

<r 2 = 0, ai > 0, 

and additionally 



73 



*3.2 



/3 2 +C (|^ 2 | 



(3.55) 



'3.3 



Moreover the admissible size of \/3 2 \ and the "O" are uniform as ui, u m and u r belong to compact sets 
of ui- and uj + . 




Figure 4: A 3-compression wave acting as a cancellation wave 

Proof of Proposition 3.6. The proof is roughly the same as for Proposition 3.4. Here we consider 

G : (/3 2 ,7s) € (-e,ef .-> ((3 2 ,j: 2 (n 3 (- l3 ,u l ),T 2 (l3 2 ,T 3 (a 3 ,u l ))). 
Again G is of class C 2 and its differential at (0, 0) is given by 

Hence we get as previously the existence of 73 cancelling a 2 and satisfying (3.55) as a consequence of 
the inverse mapping theorem. Let us now focus on the signs of o\ and 73. For 73 we have (3.55) and for 
<7i, the first-order Taylor expansion of Y,i(TZ 3 (— r y 3 ,ui),T 2 (/3 2 ,T 3 (a 3 ,ui)) gives: 



0-1=4202 +4,373 + 0(102 I') 



1 _ «|2 1 
*3,2 a | 3 a 3,3 



02+O(|0 2 | 2 ). 



(3.56) 
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Using the symmetry x < — > —x, we know from Propositions 3.2 and 3.3 that a\ 2 > 0, a\ 3 < 0, a\ 2 > 
and a| 3 < 0. Hence we can conclude that 73 < (if \(3 2 \ is small enough), but the two coefficients in 
the right hand side of (3.56) are of different signs. To conclude, we use that 63 is not too large (that 
is, we choose k small). Using Remarks 3.5 and 3.7 and adapting them in the horizontally symmetric 
situation, we see that a\ 2 = @(o"3): 033 = 0(— erf), a§ 3 = 0(— of ) and a§ 2 — ®(1)- ft follows that in 
(3.56), the second term in the parentheses is predominant over the first one for small 03. This gives the 
conclusion. □ 

4 The construction in the Eulerian case 

With the tools exposed in Sections 2 and 3, we are now in position to give our method to construct 
front-tracking approximations leading to a relevant solution for Theorem 3. We recall that front-tracking 
approximations are piecewise constant functions on the space-time domain (here R + x (0, L)), each "piece" 
on which the function is constant being polygonal. In this section, we only describe the algorithm that 
generates these approximations; we will prove in Section 6 that these approximations converge to a 
solution of the system, which will establish Theorem 3. 

The construction has some common points with the one of [21] for the controllability of the isentropic 
(2 x 2) Euler system for compressible gas, and uses some features of Bressan's front-tracking algorithm 
[6] for the generation of solutions of hyperbolic n x n systems of conservations laws with n > 3. We will 
first suppose, using the notations of Theorem 3, that: 

Ai(u ) < and A 2 (u ) > 0. (4.1) 

We will explain in Paragraph 6.1.5 how the other cases can be treated. 

As in [21], the construction consists of two successive steps. We describe these steps in separate 
subsections. A main point in the construction here is to let a strong 2-discontinuity enter the domain 
from the left side x = and to use this strong discontinuity to "eliminate" the waves inside the domain, 
using cancellation effects. As we will see, this discontinuity eventually leaves the domain through the 
right side x = L. 

We let v > a small parameter; we construct a front-tracking approximation for each such v and we 
will let v go to 0. We also let g > another positive parameter intended to go to (depending on v). 

4.1 The strong 2-discontinuity 

< 
We consider Vq such that (vq ,Uq) is an increasing (in terms of r) 2-contact discontinuity J: 

u =T 2 (a 2) Vo), <f 2 <0. (4.2) 

We require that it satisfies 

A1K7) < -Aifo), (4-3) 

which is clearly the case when \a 2 \ is small enough. Note that the velocity s of this discontinuity satisfies 

s = X 2 (u )>0. (4.4) 

This is the reference discontinuity on which the construction is based. Given such a 2-discontinuity, we 
will determine e > such that, if (1.12) is satisfied, the following construction is valid. This will allow 
us to get (1.14) by ultimately taking this reference discontinuity small. 

Now, given such a discontinuity, the approximations that we are about to construct will take values 
in the domain: 

V = B(vg;r)UB(u ;r), (4.5) 

where r > is small enough. In particular r is chosen in order that: 
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• B(vq] r) n B(uq: r) = (to simplify the discussion), 

• Pel! (in particular the vacuum is avoided), 

• any two states in B{v^\r) or in B(uq;t) determine a Riemann problem having a solution which 
avoids the vacuum, and the same is true for any "swapped" Riemann problem as defined in Sub- 
section 2.2. 

• interactions of two simple waves in B(v^ \ r) or in B(uq] r) conserve the sign in the sense of Corol- 
lary 2.1, for any permutation of the Riemann problem, and satisfy Lemma 3.1, 

• B{v^\r) C cj_, B(uQ\r) C uj+ where w_ and w+ are small enough in order for Propositions 3.1 
and 3.4 and Remark 3.2 to apply, 

• any simple wave leading a state of B(vq\t) to a state of B(u ]r) is an increasing (in terms of r) 
2-contact discontinuity with strength a 2 and speed s satisfying 

\a 2 \/2 < \a 2 \ < 2|ct 2 | and s > A 2 (u )/2, (4.6) 

• for any u in B(v^ ; r) 

Ai(tt)<Ai(«b)/2<0. (4.7) 

We will in particular choose e > in order that (1.12) implies that uq has values in B(uq; r), but e may 
have to be chosen smaller in the sequel. 

We consider A a positive number such that 

A >max |A 3 (u)|. (4.8) 

We now proceed to the construction of front-tracking approximations u v of a solution to the con- 
trollability problem; these approximations are in a first time constructed only "under/on the right" the 
strong discontinuity (in the (t,x) domain). In a second time, we resume the construction above/on the 
left of this discontinuity. 

4.2 Part 1: Construction of the approximation under/on the right of the 
strong discontinuity 

In this subsection, we describe the first part of the algorithm, which allows to construct the part of 
the approximation u" situated under/on the right of the strong discontinuity, as well as the value of u v 
immediately on the left of the discontinuity and the location of the discontinuity itself. Specifically, we 
construct the function X(t) which represents the location of the strong discontinuity in (0,T) at time 
t, and which is defined in some interval [0, Ti], T\ being the exit time of the strong discontinuity. This 
location X(t) will be an increasing function of time, depending of course on v\ to lighten the notation 
we do not make this dependence explicit. In the same time we construct the piecewise constant function 
u v on {(t,x) G [0, Ti] X [0, L] \ x > X(t)}. Moreover, we also construct the state u"(t,X(t)~) on the left 
of the discontinuity, which will be exploited in Part 2. 

In this part of the algorithm, we will suppose that all the states at points where x > X(t) belong to 
B(uQ,r) and the states u v {t,X(t)~) belong to B(vQ,r). Our convention is that the algorithm stops at 
a time when this condition starts to fail. We will prove later that the algorithm does not stop provided 
that e is small enough. 

Step 1. Approximation of the initial data and initiation of the algorithm. 

We introduce a sequence of piecewise constant approximations of the initial state (mq) in BV(0, L), 
with values in B(uo,r) and satisfying: 

TV(v,q) < TV(u ), \\uq -uolloo < \\uo- uo||oo and \\uq -u \\l 1 (o,l) < »■ (4-9) 
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Now, the algorithm to construct the approximation u v defined in M. + x (0,L) works as follows. 

a. At discontinuity point x of Uq in (0,L), we approximate the solution of the Ricmann problem 
(uq(x~),Uq(x + )) by using the accurate Riemann solver, that is by defining u u around the point (0,x) as 
the solution of the Riemann problem, where the rarefaction waves (for families 1 or 3) are replaced by 
rarefaction fans with accuracy v (described below). On the contrary, shock waves and contact disconti- 
nuities are left unchanged. 

The rarefaction fans are defined as follows: given a rarefaction wave between w_ and u + = Ti(o~i, u_), 
<jj > 0, i = 1,3, we introduce the intermediate states 



Wfe := T t [ -o-i, u_ I for k — 0, 



The rarefaction wave (2.19) is then replaced with the rarefaction fan 

{w_ for (a; — x)/t < s(u_, Wi), 
uj k for (x-x)/t e (a(wfc_i,Wfc),s(wfc,Wfe+i)), k= 1, ...,n- 1, 
m + for (a; — a7)/£ > s(w n _i, u+), 

where the shock speed for rarefactions was defined in (3.9). In other words, the rarefaction fan is 
composed of the constant states W&, separated by straight lines at shock speed s(w&, uiu+x)- 




Figure 5: The accurate Riemann solver 



b. At the point x = 0, we solve the Riemann problem (vq , Mq(0 + )); we conserve only the 2-wave and 
the 3-wave, and replace this 3-wave by a rarefaction fan with accuracy v if this wave is a rarefaction. 
The 2-contact discontinuity determines the curve X(t) for small times. At the point x = L, we consider 
that the approximation is continued with Uq(L), so that there is no Riemann problem to solve. 

After these operations, we have in (0, L) and for small times a piecewise constant function u v , where the 
constants are separated by straight lines that we call fronts; the 2-contact discontinuity originated from 
is called strong front, the other fronts being called weak. More precisely, we declare a front strong 
when it connects a state from B(vq; r) to a state in B(uq; r); it is weak when both states belong to the 
same connected component of T>. 

Moreover, all the fronts generated at this step will be called physical fronts, as opposed to artificial 
fronts which will be introduced in the next step. For each physical front separating the left state ui from 
the right state u r , there exist i E {1,2,3} (the family of the front) and oi 6 K (its strength) such that 
u r = Ti(ai,ui). 

Step 2. Extension of the solution and interactions. 

To define the approximation u v for larger t > 0, we have to explain how to extend it over points where 
two fronts meet, which are called interaction points. We do not extend any front outside of the space 
domain (0, L), so we do not have to give special rules when a front hits the boundary. 

At an interaction point (t,x), a front on the left separating the leftmost state ui from the middle 
state u m meets a front separating u m from the rightmost state u r . Of course the left front travels faster 
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than the right one. When both fronts are physical, one can write: 

u m =Ti{(ji,Ui) and u r = Tj (a'j ,u m ). (4-10) 

Remark 4.1. As in [7], we can change a little bit the speed of a front (by an amount of v at most), 
in order to avoid interaction points with more than two incoming fronts involved. Even, we can ensure 
that all times of interaction are distinct (not that this is essential). But doing so, we choose not to 
modify the speed of contact discontinuities of the second family. This is always possible since two contact 
discontinuities travelling at shock speed cannot meet, because two contact discontinuities which are not 
separated by other waves travel at the exact same speed. Also, since this can be done with an arbitrarily 
small change of speed, we avoid systematically the meeting of two rarefaction fronts of the same family 
(such a meeting does not occur naturally due to the genuine nonlinearity). Finally, we will not change 
the speed of artificial fronts (which do not meet either). 

According to the situation, the front-tracking approximation u v is extended for t > t as follows. 

A. The strong discontinuity is not involved. We suppose that none of the two fronts involved is 
the strong one. In this situation, we follow [7] (with a non-essential variant for the simplified solver). 
There are subcases. 

• Interaction with large amplitude. We suppose that both fronts are physical and that 

Wl>e. (4.ii) 

In that case we extend u v for t > t by using the accurate Riemann solver with accuracy v for 
(ui, Urn) at the point (t, x), as in the first step. However, if one of the incoming fronts (of family k) 
is a rarefaction front, we do not split the outgoing fc-wave in pieces (even if its strength is larger 
than v), and extend it as a single front sent at shock speed. 

• Interaction with small amplitude. We suppose that both fronts are physical and that 

\o l a' ] \< Q . (4.12) 

In that case we extend u v for t > t by using the simplified Riemann solver, as described now. 

— If i =/= j, that is, the incoming fronts are of different families, then i > j because otherwise 
the two fronts would not meet. The solution of the Riemann problem is approximated by 
the succession of a j-front, a i- front and an artificial front travelling at speed A. For that, 
we consider the permutation n € S3 such that 7f(l) = j, 7r(2) = i, and set £ = (1, 1, 1). We 
let a := S' r '^(uj,ti r )> Then the approximation u v is extended by a single front separating 
ui and u m :— Tj(&j,ui) travelling at shock speed, a single front separating u m from u r := 
Ti(pi, Um) and travelling at shock speed, and finally an artificial front separating u r from u r , 
and travelling at speed A. See Figure 6(a). 

— If i = j, that is, the incoming fronts are of the same family, then at least one of the fronts 
is a shock because otherwise the two fronts would not meet. The outgoing Riemann problem 
is approximated by a i-front and an artificial front as follows. Pick a permutation tt G 5 3 
such that 7r(l) = i, set £ = (1,1,1) and let a :— T^'^(ui,u r ). Then the approximation u v 
is extended by a front separating m and u r :— Ti(&i,ui) travelling at shock speed and an 
artificial front separating u r from u r , and travelling at speed A. See Figure 6(b). 

• Artificial interaction. We suppose that one of fronts is artificial. The second one is physical since 
the algorithm will guarantee that all artificial fronts under the strong front travel at the same 
speed A. Due to (4.8), the artificial front is the left one (ui,u m ); let us describe the right front 
with u r = Tj(o~j,u m ). Then one approximates the outgoing Riemann problem by a j-front and an 
artificial front as follows (we will refer to this method as the simplified Riemann solver as well). 
Pick a permutation n G S3 such that 7r(l) = j, set £ = (1, 1, 1) and let a := Y7'^{ui, u r ). We define 
u m := Tj(&j, ui), and one extends u v by a front separating ui and u m and travelling at shock speed 
and an artificial front separating u m from u r , and travelling at speed A. See Figure 6(c). 
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(c) With an artificial front 



Figure 6: Simplified solver 

Remark 4.2. We could have used the simplified solver from [7]. The (tiny) advantage here is that the 
interaction estimates enter the same framework as for the usual interactions, that is. Proposition 2.2. 
One can also notice that this simplified solver respects the fact that the interaction of two rarefactions 
of family 3 and 1 does not generate a 2-wave, and that the interaction of two shocks of family 3 and 1 
generates a 2-wave whose strength is of third order with respect to the incoming waves (recall (3.29)/ 



B. The strong discontinuity is involved. There is only one front considered strong at each time t in 

< 

this construction, of type J, separating a state in B(v , r) on the left and a state in B(uq, r) on the right, 

whose speed satisfies (4.6); otherwise the algorithm has stopped. Moreover, since in this subsection we 
are considering the approximation on the right of the strong discontinuity, at an interaction point the 
strong front is the left one. Call again ui, u m and u r the left, middle and right states. The right front 
is necessarily of the first family. Indeed, if (u m , u r ) corresponded to a physical front of the third family 
or an artificial front, it would travel faster than the 2-front (ui,u m ). And if (u m ,u r ) was separated by a 
physical front of the second family, it would travel at the exact same speed as the 2-front (ui, u m ) (recall 
Remark 4.1.) 

Hence we let u r = Ti{a\,u m ) and discuss according to the nature of this wave. 

• Interaction with a 1-rarefaction front. In that case, we use the accurate solver as described above. 
This generates a 1-rarefaction above the strong 2-discontinuity, modifies the 2-strong discontinuity 
and generates a reflected 3-wave, which is a rarefaction wave. The natures of these waves are 
deduced from the definition of r, and Remarks 3.1 and 3.2. We extend the fronts of the 1-rarefaction 
fan only for small times for the moment. This determines a new state ui on the left of the 2-strong 
discontinuity. See Figure 7. 



ut 




Figure 7: A 1-rarefaction crossing the 2-strong discontinuity 



• Interaction with a 1-shock front. In this case, we apply Proposition 3.4. We imagine that the 
3-compression wave (ui,ui) has arrived on the left of the 2-discontinuity exactly at (i,x). The 
resulting Riemann problem (v,i,u r ) for times t > t is solved by a 2-contact discontinuity and a 
reflected/transmitted wave of the third family (but no 1-wave). We know from Remark 3.8 that 
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this 3- wave is a shock. We extend u v above the interaction point (t, x) by using the accurate solver 

for the Riemann problem (ui,u r ). Consequently the front of families 2 and 3 are sent at shock 

speed. This determines ui as the new state on the left of the 2-strong discontinuity. Note that 

< 
actually, we know from Remark 3.8 that for a not too large J, the 3-wave is a shock, but this is 

not essential. See Figure 2. 

In both cases, we let X(t) follow the 2-discontinuity We do not yet extend the fronts emerging above 
the strong 2-contact discontinuity (i.e. in the domain {x < X(t)}), but we keep record of the state 
u v (t,X(i)~) above this discontinuity; this will be used in the second part of the construction. 

Note that due to (4.6), X(t) has a positive speed and eventually leaves the domain through x = L 
at some finite time T\ > with: 

Ti < j— y 4.13 

A 2 (wo) 

The first part of the algorithm ends here. 

4.3 Part 2: Construction above/on the left of the strong discontinuity 

At the end of Part 1, assuming that the algorithm is well- functioning (in the sense that it does not stop 
before T\ and generates a finite number of fronts and interaction points), we have a front-tracking ap- 
proximation defined under/on the right of the strong 2-contact discontinuity X. Let us now explain how 
we extend this front-tracking approximation u v above/on the left of the strong 2-contact discontinuity. 

Step 1. Fronts emerging from the strong discontinuity. In the construction we have left above 
the strong 2-contact discontinuity germs of 1-rarefaction waves that we would like to extend forward in 
time and germs of 3-compression waves that we would like to extend backward in time. This corresponds 
to the two following situations. At an interaction point (t, X(t)) with the strong discontinuity in Part 1, 
the state u"(t,X(t)~) on the left of the discontinuity has changed, and u_ := u"(t,X(t )~) and u + := 
u u (t,X(t )~) are connected: 

- either by u + = Ti(ai,U-) for some a\ > 0, when the incoming front on the right of X was of 
type R, 

- or by u_ = TZ^{a^, u + ) for some 03 < 0, when the incoming front on the right of X was of type S . 
See Figure 8(a). We let fronts emerge front the strong discontinuity as follows. 

• In the first situation, the rarefaction wave (ii—,u+) is treated via the usual accurate solver and 
consequently sent forward in time as a rarefaction fan. We could avoid to split these rarefaction 
waves generated by the meeting of the strong discontinuity with a rarefaction front from Part 1; 
but this has no importance. 

• In the second situation, the compression wave (u+ , U-) is split as a compression fan (with accuracy 
v) and sent backward in time. This is the equivalent for t < t of what does the accurate solver 



for rarefaction waves. To be more precise, call n := 



1 0-3 1 



and introduce the intermediate states 



u>k := 72-3(— kas/n, u_), k = 0, . . . ,n, and the propagation speeds sq := X(t ) and Si :— s(u>i-i,uii) 
for i = 1, . . . , n. Then u v is set on the left of X locally at (I, x) as 

u u (t,x) — u>i for t < t, x<X(t) and = — € [s,, Sj+i], i = 0, . . . ,n — 1, 

u v (t, x) — u + for t < t, x < X(t) and = — > s n or for x < X(t) and t >t. 

It is clear that S\ > \^{u^) > Aa(u-) = sq and that Sj+i > Si for % = 1, . . . , n— 1 since compression 
waves satisfy Lax's inequalities: 

X i (TZ i (a l ,u)) < s(u,TZ t (ai,u)) < A.,(u), for i = 1,3, Oi < 0. (4-14) 
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x = 




x = L 




t = t = Ti 

(b) Beginning of Part 2 



(a) At the end of Part 1 

Figure 8: Junction of the two parts of the construction 



Step 2. Extension of the fronts. Now we have to explain how we extend these fronts and complete 
the approximation u v on the whole M + x (0, L). The main point is to use L — x as time variable; we are 
led to an initial-boundary value situation with a moving boundary, see Figure 8(b). In order to avoid 
the confusion with the actual time variable t, we will systematically refer to this variable L — x as the 
pseudo-time. The "initial data" on {L} x [Ti,+oo) is u v (T^ ,L~) (this state on the left of the strong 
discontinuity was determined during the first part of the construction), and the entering waves from the 
moving boundary {(t,X(t)), t £ [0,Ti]} are the germs mentioned above. 

Now we start from the state m 1/ (T 1 + , L~) at the pseudo-time L — x = 0, and let L — x increase. When 
an interaction point with the strong discontinuity obtained in Part 1 is met, we let the fronts enter the 
domain as described in Step 1. Note that all these fronts evolve forward according to the pseudo-time. 
We have to explain how to extend the approximation u v when an interaction point inside the domain 
x < X(t) is met. As we will see, only one case can actually occur. 

• Interaction of a 1-raref action front with a backward 3- compression front. Assume as in Figure 9 that 
at some pseudo-time L — x = L — x and actual time t = t, a 1-rarefaction front u m = Ti(<7i, ui), o~± > 0, 
meets a 3-compression front u m = TZ^(a3, u r ), 03 < 0. Then one solves the swapped Riemann problem 
(see Subsection 2.2): 

u r = TZi(a[, •) o K 2 {o-' 2 , ■) o ^3(^3, •) m. 

Using Lemma 3.1, we see that a' 2 = 0. The fact that the waves of families 1 and 3 conserve their nature 

(C and R) across the interaction point, or in other words that a[ > and <j 3 < 0, is a consequence of 
Corollary 2.1 and the definition of r. We denote 

U m = K3(o-3, Ul)- 

The approximation u v is extended for further pseudo-times as a backward 3-compression wave separating 
u\ and u m travelling at shock speed and a forward 1-rarefaction front separating u m and u r and travelling 
at shock speed. 

Precisely the wave pattern for L — x > L — x is locally given as follows: 



u v (t,x) 



ui for £ < — t= r, 

u m for t^ — r < £ < 

u r for £ > l 



s(u m ,u r ) ' 



with £ = 



t-t 



We will refer to this construction as the side simplified solver. 
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L — x 



Figure 9: Simplified solver for side interactions 



The important fact here is that there is no other interaction occurring in this part, other than of the 
type described above. Let us justify this. 

• In the algorithm above, only forward 1-rarefaction fronts and backward 3-compression fronts enter 
from the boundary. Since there are no front initially (for time L — x = 0) and since only 1- 

rarefaction fronts and backward 3-compression fronts emerge from a R/C interaction point, there 
are only forward 1-rarefaction fronts and backward 3-compression fronts in the domain as long as 
no interaction of another type occurs. 

• In the algorithm above, no front goes back to the strong 2-discontinuity: the forward 1-rarefaction 
fronts because they go forward in time at negative speed, the backward 3-compression fronts 
because, in the usual sense of time, they travel faster than the 2-discontinuity. 

• There are no interactions of fronts within a family since 1-rarefaction fronts traveling forward are 
going away one from another, so do 3-compression fronts when going backward in time (see (4.14)). 

Hence the description of the algorithm is complete. 

Remark 4.3. The advantage of using 3-compression waves is that their interactions with 1 -rarefactions 
do not generate a wave in the second family. However we could have used fans of small 3-shocks to 
replace the 3-compression waves fan. The (small) cost would have been the appearance of artificial fronts 

(travelling to the left) at each R/S interaction point. But since in that case, the resulting artificial front 
is of third order with respect to the sizes of the incoming waves, the estimate of the total strength of these 
artificial fronts would have been relatively easy. 

5 The construction in the Lagrangian case 

The construction in the Lagrangian case, still relying on a front-tracking algorithm, is of different nature 
than in the Eulerian case, since, obviously, one cannot make a 2-contact discontinuity travel through 
the domain. Here we will use two successive strong shocks: first, a 1-shock crossing the domain from 
right to left, and then a 3-shock crossing the domain from left to right. There will be three parts in the 
construction: first under/on the left of the strong 1-shock, then above/on the right of the strong 1-shock 
but before the entrance of the 3-shock, and finally after the entrance of the strong 3-shock. A main 
difficulty here is to eliminate the 2-contact discontinuities, since they have zero characteristic speed and 
hence do not propagate to the boundary. 

We recall that the construction below is also valid in the Eulerian case when 7 < f , with minimal 
changes. 

As before we let v > a small parameter and construct a front-tracking approximation u v for each 
v small. We also let g > another positive parameter to be determined (depending on v). 
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5.1 The two strong shocks 

We consider v^ such that (uo,Vq) is 1-shock S: 

v£ = Ti(ai,u ), cti < 0. (5.1) 

Its velocity si satisfies 

S!<Ai(u )<0. (5.2) 

Then we consider uj~ such that (v^jVq) is 3-shock S: 

Vq =T 3 (a 3 ,v^), 0=3 < 0. (5.3) 

We suppose that a 3 is small enough for Proposition 3.6 to apply. The velocity S3 of this shock satisfies 

s 3 > A 3 «) > 0. (5.4) 

Given these shocks, we introduce the domain: 

V = B(u ; r) U B(v+;r) U B(v^;r), (5.5) 

and choose r > small enough in order that: 

• B(uo;r), B(vQ\r) and B(v^;r) are disjoint, 

• Pel! (in particular the vacuum is avoided), 

• any two states in B(uo;r), in B(vq ; r) or in B(uj~;r"), determine a Riemann problem having a 
solution which avoids the vacuum, and the same is true for any "swapped" Riemann problem as 
defined in Subsection 2.2. 

• interactions of two simple waves in B(uo;r), in B(vq ;r) or in _B(u^;r), conserve the sign in the 
sense of Corollary 2.1, for any permutation of the Riemann problem, and Remarks 2.2 and 2.3 
apply, 

• Propositions 3.3 and 3.5 apply with B(u ;r) C u;_ and B(vq]t) c uj + ; Proposition 3.6 applies 
with B(vq ;r) C w_ and B(v^;r) C cj + , 

• any simple wave leading a state of B(uo;r) to a state of B{vq\t) (resp. a state of B(v^;r) to a 
state of B(vq ;r)) is a 1-shock (resp. a 3-shock) with strength <j\ (resp. 173) and speed s satisfying 

|ffi|/2 < I o-i I < 2|cti| and s < Ai(u )/2, (5.6) 

(resp. |a 3 |/2< M < 2|a 3 | and s > A 3 «)/2,) (5.7) 

and moreover Proposition 3.6 applies to any such 3-shock, 

• for any u in B(vq ;r) 



and for u in B(v 1 ; r) 
We introduce A satisfying: 



A3(«) > A 3 «)/2 > 0, (5.8) 

Ai(«)<Ai(«r)/2<0. (5.9) 



A > max |Ai(u)|. (5.10) 
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5.2 Part 1: Construction below/on the left of the strong 1-shock 

In this first part of the construction, we describe the design of u v under/on the left of a strong 1-shock, 
which enters the domain from x — L at time 0, and eventually leaves the domain through x — 0. Together 
with u v we construct the location of this shock which is described by the function X\ (t) , defined in some 
interval [0, Ti], T\ being the exit time of the strong 1-shock. Thus the piecewise constant function u v 
is determined on {(£, x) <G [0, Ti] x [0,L] | x < X\(t)} after this part. This part of the algorithm also 
provides the value u u (t,X^~(t)) immediately on the right of the strong 1-shock. 

Again the algorithm is supposed to generate states belonging to B(u , r) in {(£, x) £ [0, Ti] X [0, L] | x < 
Xi{t)} and states belonging to B(vq ,r) on (t,X^(t)); we consider that it stops as soon as it should 
generate another state. We will prove later that provided that e is small enough, the algorithm does not 
stop. 

Step 1. Approximation of the initial data and initiation of the algorithm. 

As in Section 4, we initiate the algorithm by introducing a sequence of piecewise constant approxi- 
mations of the initial state (uq) in BV(0, L), with values in B(uq, r) and satisfying (4.9). Then we start 
the algorithm as follows. 

a. At a discontinuity point x of Uq in (0, L), we approximate the solution of the Riemann problem 
(uq(x~), Uq(x + )) by using the accurate Riemann solver, exactly as in Section 4. 

b. On the right point x = L, we consider the Riemann problem (uq(L~),Vq ) and approximate its 
solution by using the accurate Riemann solver, and conserve only the 1-wave, which is a 1-shock due to 
the restrictions on T>. On the contrary, x = is considered a continuity point. 

This first step determines the various states of u v and the location X\(t) of the 1-shock for small times. 
As before, in order to define u v for later times, one must describe what happens at interaction points. 
As in Section 4, we do not extend any front outside of the space domain (0, L), so we do not give rules 
concerning a front hitting the boundary. 

Step 2. Extension of the approximation and interactions. 

At an interaction point (t, x), a front separating the leftmost state ui from the middle state u m meets 
a front separating u m from the rightmost state u r , which we write again (4.10) when both fronts are 
physical. Of course, the left front has a larger speed than the right one. Again we can change a little 
bit the speed of a front (of an amount at most of v) , in order to avoid interaction points with more than 
two incoming fronts involved. We do not modify the speed of the contact discontinuities (of family 2), 
though, and we avoid the meeting of rarefaction fronts of the same family. 

According to the situation, the front-tracking approximation u v is extended over an interaction point 
as follows. 

A. The strong 1-shock is not involved. In this case, we use the exact same strategy as in Section 4 
that is: 

• Interaction with large amplitude. If both incoming fronts are physical and (4.11) is satisfied, then 
we extend u v by using the accurate Riemann solver with accuracy v. 

• Interaction with small amplitude. If both incoming fronts are physical but (4.12) is satisfied, we 
extend u" by using the simplified Riemann solver as in Subsection 4.2. However here we set —A 
as the speed of artificial fronts, which are therefore placed on the left of the outgoing fronts. This 
means the following: 

— if the two incoming fronts are of different families i > j, we use the permutation 7r <G 5.3 
determined by tt(2) = j and 7r(3) = i, £ = (1, 1, 1), we consider the corresponding swapped 
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Ricmann problem, and extend the approximation via an artificial front at speed —A, a j-front 
and a j-front at shock speed, 

— if the two incoming fronts are of the same family i — j, we use a permutation n £ S3 such 
that 7r(3) = i, £, = (1, 1, 1), consider the corresponding swapped Riemann problem and extend 
the approximation via an artificial front at speed —A and a z-front at shock speed. 

Hence the situation would be described by Figure 6 after a vertical symmetry. 

• Artificial interaction. If one of the fronts is artificial, we use the simplified Riemann solver, as we 
have described it in Section 4, taking the speed —A of artificial fronts into account. Again one can 
think of Figure 6 after a vertical symmetry. This amounts to considering a permutation n £ S3 
such that 7r(3) = j, where j is the family of the incoming physical front, and £ = (1,1,1). 

B. The strong 1-shock is involved. This is where the strategy changes. We consider the interaction 
of the weak physical front (ui, u m ), let us say, 

u m =Ti(ai,ui), (5-11) 

with the strong 1-shock (u m , u r ) (which is the continuation of the 1-shock emerging from x — L at initial 
time). The weak wave is on the left since we construct the approximation u v under/on the left of the 
strong 1-shock. There are no interactions between the strong 1-shock and artificial fronts, since the latter 
are faster. 

There are six cases according to the family i of the weak wave and to its nature (<Xj > or <7j < 0) . 
These six cases are gathered in two groups. 

*- > -> 

• Group I: S , J and R. In this group, the incoming weak front is either a 1-shock, a decreasing 

2-contact discontinuity or a 3-rarefaction. In that case, we use the usual accurate Riemann solver 

for the outgoing waves. We do not yet extend the outgoing fronts further in time, except for the 

1-shock which describes X\. According to Proposition 3.3 and the definition of r, the outgoing 

> '- 

wave in the family 2 is a J contact discontinuity and the outgoing wave in the family 3 is a R 

rarefaction. (This is the reason which brings together these incoming waves in this group.) This 

determines a new state on the right of X\ . 

• Group II: R, J and S . In this group, the incoming weak front is either a 1-rarefaction, a in- 
creasing 2-contact discontinuity or a 3-shock. Here we use a correction wave. Indeed, we apply 
Proposition 3.5, and imagine that a 1-compression wave (u r , u r ), determined as in this proposition, 
has arrived at (t, x) on the right of the strong 1-shock and solve the resulting Riemann problem 
(ui, u r ). We do not extend the outgoing fronts yet, except for the outgoing 1-shock which continues 
X±. Taking this additional 1-compression wave into account, the outgoing wave in the family 2 is 

a J contact discontinuity and the outgoing wave in the family 3 is a R rarefaction as well, as a 

consequence of Proposition 3.5 and the definition of r. If this additional 1-compression wave was 

< ■<- 

not there, we would obtain a J contact discontinuity and a S shock in families 2 and 3 respectively. 

This determines the new state on the right of X\. 

These two situations allow in particular to extend the strong 1-shock and the function X\(t) further in 
time, and to keep track of the state on the right of the strong shock (as long as the algorithm has not 
stopped). 

Now, if the algorithm has generated only a finite number of interactions and if it has not stopped, 
then, due to the definition of r, the strong 1-shock leaves the domain through x = at some time T\ 
satisfying 

T x < jy^m- (5-12) 

We represent the situation at the end of Part 1 in Figure 10. 
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Figure 10: End of Part 1 

5.3 Part 2: Construction between the two strong shocks 

At the end of the first part of the construction and assuming that the algorithm is well-functioning, we 
have obtained a front-tracking approximation u v defined under the strong 1-shock. Let us now explain 
how we extend u v above/on the right this shock. 

To complete the approximation u" we extend the 2-discontinuities and the 3-rarefaction fronts forward 
in time and the additional 1-comprcssion waves backward in time. For that purpose, we will not quite 
use a; as a new time variable (though this gives the main idea), but rather the variable 



t) 



it. 



with i > chosen small enough so that i is smaller and strictly separated from {|Ai(w)|, ^(u), u £ T>}. 
In particular the lines $ = constant have a slope strictly separated from the one of characteristics or 
shocks of the three families as the states belongs to T>. For this new time variable, all the fronts emerging 
from Xi (including the 2-contact discontinuities) go "forward" . As in Section 4, we will refer to d as the 
pseudo-time to avoid confusion with the actual time t. 

We represent the situation at the beginning of Part 2 in Figure 11. 

The algorithm here is supposed to generate states in B(vq ,r) on {x > Xi(t)}; as in the first part, 
we consider that it stops as soon as it generates another state. 

Step 1. Fronts emerging from the strong 1-shock. In the construction we have left above/on the 
right of the strong 1-shock germs of 1-rarefaction and 2-contact discontinuities to be extended forward 
in time and germs of 3-shock fans to be extended backward in time. See Figure 10. More precisely, at 
a point from where the fronts emerge, let us say (t,x) = (t,Xi(t)), the two states u_ := u v (t ,x~) and 
u + := u" '(t ,X~) are connected through 

u- =1li(-<ri,-) °T 3 (a 3 ,-) oT 2 (a 2 ,-)u + , 

with <7i < (potentially there is no additional 1-compression wave), 03 > and a 2 < 0- We introduce 
the intermediate states u^ := T%{a<x-, ') u + anc ^ u ™ := ^M ^ ') u m- 
We let fronts emerge from the strong discontinuity as follows. 



The contact discontinuity (u + ,ti^) and the rarefaction wave (u^ 
accurate solver and sent forward in time. 



are treated via the usual 



• The compression wave (u_,w^J is as in Section 4 split in fronts of size at most v which are sent 
backward at shock speed. Again, from Lax's inequality (4.14), these fronts emerge indeed from the 
strong 1-shock inside the zone {x > X(t)}. 

Step 2. Extension of the fronts. Now, the rest of the algorithm consists in extending fronts across 
interaction points; we do not extend fronts outside of the space domain (0, L). In order to construct u v 
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Figure 11: Beginning of Part 2 

above/on the right of the curve X±, we start from u v {T\, + ) for (t, x) £ (T\, +00) x {0} (this state was 
determined in the first part of the algorithm) and progress with the pseudo-time variable ■&. When an 
emergence point on X\ is met, we extend the fronts outgoing from X\ as described before (these fronts 
all go forward according to the time variable •&). We have to explain what we do at interaction points 
inside {x > Xi(t)}. 

First, as will be clear from the algorithm, only fronts of the following nature will be produced: 3- 
<- > 

rarefaction fronts R (going forward in time t), 2-contact discontinuities J (going forward in time t) and 

1-compression fronts C (going backward in time i). It follows that there will be no interaction between 
these fronts and the strong 1-shock generated by the first part of the algorithm: the waves R and J 

because they go forward at a non-negative speed, the waves C because they go backward in time and 
satisfy Lax's inequalities. It follows also that there will be no interactions between fronts of the same 

family: the waves R because they go forward and spread (because they have positive strength and 

because of (2.7)), J because they all go forward at the exact same speed (that is zero), the waves C 
because they go backward and satisfy Lax's inequalities. 

The extension of u v beyond an interaction point depends on the nature of the incoming fronts, which 
are all weak waves. 

• Interaction of R and C. We consider the situation where a backward 1-compression front (ui,u m ) 
with ui above u m in the (t,x) plane, meets a forward 3-rarefaction front (u m ,u r ) with u m on the 
left of u r . This is described in Figure 12(a). One has ui = lZi(a, u m ), a < and u r = T$(fi, u m ), 
P>0. 

In that case, we use the same type of solver as in Subsection 4.3. Precisely we solve the swapped 
Riemann problem 

u r = TZi(-a[,-) o 72.2(02,-) K 3 ((T3,-)mi. 

Due to Lemma 3.1, one has <r 2 = 0. Moreover due to Proposition 2.2 and the definition of r, one 
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has a'i < and a' 3 > 0. Setting u m := T 3 (a' 3 ,ui), we extend the approximation for further •& via 
a single forward 3-rarefaction front (ui,u m ) and a single backward 1-compression front (u m ,u r ), 
both sent at shock speed. 

> ->■ 

• Interaction of a contact discontinuity J and a rarefaction R. We consider the case where a for- 
ward 3-rarefaction front (ui,u m ) with ui on the left of u m in the (t, x) plane, meets a forward 
2-discontinuity (u m ,u r ) with u m on the left of u r . This is described in Figure 12(b). One has 
u m = TZ 3 (a, ui), a > and u r = T 2 (j3, u m ), /3 < 0. 

In that case, we apply Corollary 2.2 to System (1.3), with k = 1. In other words, we seek to cancel 
the "outgoing 1-wave" (outgoing for the usual time t). We deduce the existence of 7, such that if 
one sets u r := 1Zi{^,u r ), then Tt\{ui,u r ) = 0, and 

|E2(uj,Ur) -0\ + |S 3 (ui,«r) - a| = 0(\a\\p\), 7 = -otfi li ■ [r 3 , r a ] + 0((|a| + |/3|)|( 



Now 



?i-[r 3 ,r 2 ] = -4~- ~ n. 



9r 4r 

Hence we deduce that 7 < 0, that is, (u r ,u r ) is a 1-compression wave. Due to Corollary 2.1 and 

> -> 
the definition of r, the other two resulting waves are of type J and R. 

We extend the approximation u v over the interaction point, for further $, by the four states u/, 

> -* *~ 

iii, Ur and u r separated by the three (single) outgoing fronts J, R and C traveling at shock speed; 

> -1, . i_ 

this is forward in time for J and R, and backward in time for C. 




1? = est 

(a) Interaction R/C 




J § = est 
(b) Interaction J/R 



6 

x 




■d = cst 



> '- 

(c) Interaction J/C 



Figure 12: Side interactions 



Interaction of a contact discontinuity J and a compression wave C . We consider the case where 
a backward 1-compression front (u m ,ui) with u m below ui in the (t,x) plane, meets a forward 
2-discontinuity [u m ,u r ) with u m on the left of u r . This is described in Figure 12(c). One has 
u m = Tli(-a,U[), a < and u r = T 2 (/3,u m ), (3 < 0. 

In that case, we apply Corollary 2.3 to System (1.3), with k = 3. In other words, we seek to cancel 
the "incoming 3-wave" . We deduce the existence of a, such that 



|<7i - a\ 



u r = Ki(-ai, •) o T 3 (er 3 , •) o T 2 (a 2 , ■) Ui, 
\a 2 -fi\ = 0{\a\\p\), a 3 = a{3£ 3 - [n,r 2 ] +0((\a\ + 



Here 



or 4r 



«3 • [ri,r 2 \ ■■ 

Define u m := T 2 (o~ 2 ,u{) and u r = T 3 {o~ 3 ,u m ) = TZi(o~i,u r ). Using the definition of r, we deduce 

> 
that <7i < 0, o~ 2 < and a 3 > 0, that is, (u r ,u r ) is a 1-compression wave, (ui,u m ) is a J wave and 
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(u m ,u r ) is a 3-rarefaction wave. We extend the approximation u" over the interaction point, for 

further i?, by the states ui, u m , u r and u r separated by the three (single) outgoing fronts J, R and 

C traveling at shock speed, the fronts J and R moving forward in time, and the front C backward 
in time. 

The description of the algorithm for this second part is complete since no fronts are created other than 
C, J and R, and the possible interactions between all these types of fronts were covered. 

Now we claim the following. 

Lemma 5.1. Assuming that the algorithm generates an approximation for all times (with a finite number 

> 
of interaction points), there are only fronts of type J present in the domain (0, L) for times t > T 2 , 

T 2 :=T 1 + -^—. (5.13) 

A3W) 

Proof of Lemma 5.1. Consider above X\, at some time t , a 3-rarefaction front or a 1-compression front 
a; call x& its position at time t , and let f? 1 = x& + ti 1 the pseudo-time associated to the point (t , %). 
We have the following algorithm to get back (according to the pseudo-time) to the "origin" of this front, 
that is, an interaction point on X± where we consider the front to be coming from. 

Step 1. Given a front a, we go to the "earlier" interaction point (according to $), that is, 

- if the front is a 3-rarefaction, we go back in time t to the previous interaction point, 

- if the front is a 1-compression front, we go forward in time t to the next interaction point. 

Step 2. At the interaction point, we discuss according to the nature of this point (see again Figure 12): 

- if it is a C / R interaction point, follow for earlier •& the incoming front of the same family as a and 
go to Step 1, 

- if it is a C/J interaction point, follow for earlier $ the incoming C front and go to Step 1, (whether 
a is a C front or not), 

- 1 > -* 

- if it is a Rj J interaction point, follow for earlier $ the incoming R front and go to Step 1, (whether 

a is a R front or not), 

- if it is an interaction point with the strong shock on X\, stop here. 

Note that this algorithm stops, since we go from an interaction point to another with decreasing •d, and 
there is no front coming from x = 0, t > T-y. 

Now following the front a from its origin on X\, say (t°, Xi(t )), to (t 1 , Xa) in increasing pseudo-time, 
there are pseudo-time-intervals [$2i)$2i+i] where a is a 3-rarefaction front going forward in time t, and 
pseudo-time-intervals [#2i+i ^21+2] where a is a 1-shock front going backward in time t. Call (U,Xi) the 
interactions points corresponding to pseudo-times $j. 

The real time t increases during the pseudo-time-intervals [$2i,$2i+i\i and decreases during the 
pseudo-time-intervals [^21+1,^21+2]; the position x progresses for all pseudo-time-intervals \di,di + i\. 
Hence, recalling that the 3-characteristic speed is bounded from below by A3(i*o")/2, we deduce: 

t -t < > hi+i - hi < . , + , / , a?2i+i - %2i < - , xr i 

which proves the claim. D 

We consider the second part of the algorithm to stop at time T2, where the third and last part of the 
algorithm begins. 
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5.4 Part 3: Construction using the strong 3-shock 

Let us describe how we extend u u for times t > T 2 . The idea here is to let a 3-strong shock based on the 

reference shock (v± , Vq) enter the domain from x = and eventually leave the domain through x = L. 

This 3-strong shock will allow us to get rid of the remaining fronts, which, according to Lemma 5.1, are 

> 

J fronts. 

The algorithm here is supposed to generate states in B(vq ,r) below/on the right of the 3-strong 

shock and states in B(v± , r) above/on the left of the 3-strong shock; it stops otherwise. 

Step 1. Under/on the right the strong 3-shock. 

> 
Initial data. At time T 2 , we have vertical J fronts in (0, L) coming from the second part of the 

algorithm. To these fronts we add at x — the solution of the Riemann problem (v± , u u (T2, + )), from 

which we conserve only the 3-wave. Using the definition of r, this 3-wave is a (strong) 3-shock. We call 

A 3 (£) its position at time t. Now the goal in this step is to construct u v for t > T 2 , x > X 3 (t), together 

with the position A3 of this strong 3-shock and with the state u v (t, X 3 {£)~) on the left of this shock. 

Interactions. Since J fronts do no interact between themselves, having all zero speed, we only need to 

specify what happens at an interaction point between the strong shock S and a weak front J. 

> 
We suppose that the strong 3-shock (ui, u rn ) meets a decreasing 2-contact discontinuity J with states 

(u m ,u r ). We apply Proposition 3.6. We deduce that there exists 73 < such that, assuming that a 

3-compression wave {ui,u{) with strength 73 arrives exactly at this interaction point from the left side 

of A3, there is no outgoing J wave, and the outgoing wave for the first family is a rarefaction one. We 

extend A3 by following the outgoing 3-shock at shock speed, and consider the state on the left of the 

outgoing 3-wave in Proposition 3.6 as the state on the left of A3. See Figure 4. 

This allows to construct the approximation on the right of A3. Since, assuming that the algorithm has 

not stopped, all the states on the right (respectively left) side of A 3 belong to B(vq , r) (resp. B(v^ , r)), 

the speed of the shock satisfies A 3 > X 3 (vq)/2, so the strong shock leaves the domain at some time T 3 

satisfying: 

T 3 <T 2 + -^—. (5.14) 

Step 2. Above/on the left of the strong 3-shock. 

We have left to extend u u above A3. For this, we follow the same method as in Section 4, see Figure 8. 
The only difference here is that the strong discontinuity is no longer a 2-contact discontinuity but a 3- 
shock; but this does not intervene since there are no new interactions with the strong discontinuity. 

Hence, we use L — x as a pseudo-time variable. We start from the state u v (T 3l L~) at pseudo-time 
L — x — (on the "space domain" which is originally [T 3 , +00)). Wc let the pseudo-time L — x progress 
until one meets an interaction point inside the domain, or on the boundary (L — x, t) = (L — X 3 (t), t). 

Fronts emerging from the boundary. At an interaction point on A3 coming from Step 1, we have 
U[ = u u (t~ , A 3 (t)) and u r ~ u"(t + . X 3 (t)) connected via a 3-compression wave and a 1-rarefaction wave: 

u r = Ti(cri,M m ), ui = ^3(73, u m ), cti > and 73 < 0. 

We extend the approximation u u for larger L— x by tracing between u m and ui a backward 3-compression 
fan with accuracy v replacing the actual 3-compression wave; as before we let the fronts go at shock speed 
backward in time t, that is forward with respect to L — x. We approximate the rarefaction wave (u m , u r ) 
by using the accurate solver (splitting it in pieces no larger than v) and sending the corresponding fronts 
forward in time at shock speed. 

Interactions. Interactions of two fronts C and R are treated using the side simplified solver, exactly as 
in Subsection 4.3 (see in particular Figure 9); again due to Lemma 3.1 no 2-wave appears here, and due 
to Proposition 2.2 the new waves are of the same nature as the incoming ones. There are no interactions 
between waves of the same family (the rarefaction fronts go forward in time, the compression fronts go 
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backward in time and satisfy Lax's inequality) and no interaction with the strong 3-shock (rarefaction 
fronts go forward at negative speed, shock fronts go backward in time and satisfy Lax's inequality.) 

This ends the algorithm in the Lagrangian case. 

6 Proofs of the main results 

In this section, we establish Theorems 1 and 2, starting with the proofs of Theorems 3 and 4. For that, 
we prove that the wave front-tracking algorithms described above are well-functioning, in the sense that 
they generate an approximation u v defined for all times, with a finite number of fronts and interaction 
points. There could be otherwise at some time an accumulation of interaction times. At the same time, 
we prove estimates on the sequence u u which will allow to extract a converging subsequence and to 
prove that the limit is a suitable solution of the problem. Several parts are done by adapting [7] to our 
situation. 

We will denote as in [7] the fronts by greek letters. When a is a front, a a G K denotes its strength, 
k a 6 {0,1,2,3,4} its family (with the convention that 4 describes the family of artificial fronts in the 
case of Theorem 3, describes the family of artificial fronts in the case of Theorem 4). 

6.1 Eulerian case: proof of Theorem 3 

There are several successive steps. 

6.1.1 BV estimates 

We first prove uniform BV estimates on the approximation u u (as long as it is well-defined and all the 
states belong to T>, which a posteriori will be proven to be all times). This is an adaptation of Glimm's 
argument [24], but here Glimm's functionals have to be defined along curves which are suited to the 
geometry of the construction. 

We introduce for each time t two curves Tj and Tf, and for each x £ [0,L] a curve C%, all these curves 
being drawn inside R + x [0, L\. Our goal is to bound the total variation of the approximation u" along 
these curves. These curves depend on i/, but to lighten the notation we do not make this dependence 
explicit. Recall that the 2-strong discontinuity X has a positive velocity. 

• Given t > 0, we define the curve T] as the part of the curve X describing the 2-strong discontinuity 
for times in [0,£], glued with the horizontal curve {t}x [X(t), L]. One should have the representation 
that the part of X that is considered is X(t) + , that is, the right side of the discontinuity. See Figure 
13(a). We do not consider Tj for t larger than T\. 

• The curves C^. are obtained by gluing the part of the curve X describing the 2-strong discontinuity 
situated in the space interval [0, a;], with the vertical line segment [X~ 1 (x), +oo) x {x}. Here the 
portion of X considered is the one on the left side. See Figure 13(b). 

• The curves Tf are obtained by gluing the vertical line segment [t, +oo) x {0}, the horizontal line 
segment {£} x [0, X(t)] and the part of the curve X describing the 2-strong discontinuity for times 
in [t, Ti]. Again one has in mind that the portion of X is considered on the left side. See Figure 
13(c). For t > Ti, the curve rf is only composed of the vertical line segment [t, +oo) x {0} and 
the horizontal line segment {t} x [0, L]. 

We follow rj from left to right according to the variable x, we follow C\ from bottom to top according 
to the variable t and finally we follow Tf by first following the vertical line segment from top to bottom, 
then following the rest of the curve from left to right. We will say that we follow these curves from "left 
to right" when we follow them in this way. Given two points on one of these curves, this gives a meaning 
to "one being on the left of the other" . 

We consider the following Glimm functionals, i = 1,3: 

V\t)= Y, { 1 + ^ l i6 lka S Xa>x{t) )\a a \ and Q l (t) = £ KIM, (6.1) 

a front . a fi f ront cutting T\ 

cutting 1 t a a pp r oaching /3 
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(a) The curve Tj 





(b) The curve Cl 

Figure 13: Curves for the BV estimate 



(c) The curve rf 



and 



V 2 (x) 



^2 Wa\ and Q 2 ( X ) 



E 



kal \(T/}\, 



(6.2) 



a front 
cutting C 



a,/3 front cutting C x 
a approaching /3 



Let us give some precisions: 



• We consider that a front crosses T\ in its part coinciding with X only when it leads into/emerges 
from X on its right, not when it emerges from X on its left. In the same way, only fronts emerging 
from X on its left can cross rf or C% (assuming that they touch the correct part of X), not a front 
leading into/emerging from X from its right. 

• Note that a 3-front can cut rj twice (once when emerging from X(t) on its right, once cutting the 
horizontal part of rj-); it this case we count this front twice. In the same way, a 1-front can cut Tf 
(and be counted) twice. 

• We define a and j3 cutting T] as approaching, when, a (of family i) being on the left of j3 (of family 
j), one has i > j or i = j and at least one of a or /3 is a shock. Here artificial fronts are considered 
of family 4. 

• We consider a and f3 cutting C^ as approaching when, a being on the left of j3, the couple (a, f3) 
is of type (R,C). 

• We consider a and /3 cutting T^ as approaching when, a being on the left of /3 in the sense given 
above, the couple (a, /3) is of type (C, R). 

• k > 1 is a constant to be determined. It is used to put a slightly different weight when the front 
a corresponds to i = 1 (a front under the strong shock) , k a = 1 (a front of the first family) and 
x a > X(t) (the front cuts r£ on its horizontal part, so that 5n is actually redundant). This is to 
take the reflections of waves by the strong discontinuity into account. 

Now applying Glimm's method we prove the following lemma. 

Lemma 6.1. For k > and K > large enough, ifTV(uo) is small enough, then 

F 1 ^) := V 1 ^) + KQ x {t) is non-increasing. (6.3) 

Proof of Lemma 6.1. The function .F 1 is piecewise constant in time. Actually it is modified at an inter- 
action time or at a time where a front leaves the domain (0, L). In the latter case, both the functionals 
defining F 1 decrease, so that we only need to understand what happens at interaction times. There are 
two distinct types of interactions. 
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• Interaction of weak waves. When at the time of interaction t, two weak waves a and /3 interact, 
then one can perform the classical analysis (relying on Proposition 2.1 when the accurate solver is 
used, and on Proposition 2.2 when the simplified solver is used): 

V 1 ^) ^V 1 ^-) +C lK \a a \\ap\, and Q 1 ^) < Q 1 ^') - \a a \\a fi \ + C 2 \a a \\a ? \F 1 (t'). 

As a consequence, given k > 1, there exist K > and S e (0, 1) such that if F 1 ^) < 5, then one 
has F 1 (i + ) < F 1 (t _ ) when crossing an interaction of weak waves. Even, in that case, we can have 

F\t+)~F\t-)<-\a a \\a p \. (6.4) 

• Interaction with the strong wave. Let us consider a time of interaction t, when a weak front a meets 
the strong discontinuity. The front a is necessarily of the first family. Whether it is a rarefaction 
front or a shock, its interaction with the strong discontinuity will result in a reflected 3-wave /3 
which crosses Tj twice for times just after t. Moreover one has the estimate (see Propositions 3.1 
and 3.4) 

M<c 3 M. 

Since after t, the front a does no longer cut T] on its horizontal part, it follows that 

V\t + )<V\t-) + 2C 3 \a a \-n\a a \, and Q\t + ) < Q l (r) + C 4 \a a \F l (r). 

Hence one can find n > such that if F x (t _ ) < 1, then one has F x (t + ) < F 1 ^) when crossing 
an interaction of a weak wave with the strong discontinuity. Even, in that case, we can have 

FHO-^O^-KI- (6-5) 

The above analysis allows to find 5 > such that if F 1 (Q) < 5, then F 1 is non-increasing. Since one has 
F 1 ^) < C 5 (kTV(u ) + KTV(u ) 2 ), one deduces the claim. □ 

The same method applies to C 2 and T 3 , which leads to the following statement. 

Lemma 6.2. 1. If K > is large enough and ifTV(uo) is small enough, then 

F 2 (x) := V 2 (x) + KQ 2 {x) is non- decreasing, (6-6) 

F 3 (t) := V 3 (t) + KQ 3 (t) is non-increasing. (6.7) 

2. For some C > 0, one has 

F 2 {L)<CF 1 (T 1 ), (6.8) 

F 3 {0) <F 2 (0) + F 2 (L). (6.9) 

Proof of Lemma 6.2. The first part of Lemma 6.2, that is (6.6) and (6.7), is analogous to Lemma 6.1, 
relying on Proposition 2.2. Actually, it is even simpler since there are no artificial fronts above the strong 
discontinuity and no wave interact with the strong front in the second part of the construction. Hence 
we do not need a n5ik a here. Note that V 2 and Q 2 do not change when x corresponds to an interaction 
point with the strong discontinuity and only weak interaction points may affect F 2 . In the same way at 

a time where a front exits through x = 0, the functional V 3 and Q 3 do not change (for C fronts) or 

decrease (for R fronts) . Also at a time of interaction with the strong shock, the functionals V 3 and Q 3 

do not change (if there are no C fronts arriving there) or decrease (otherwise) . 

Concerning (6.8), the values F 2 (L) and i 7 ' 1 (Ti) measure the total strengths of the waves on the left 
and right sides of X, respectively (remark that no front leave the domain through (Ti,+oo) x {L}). 
To get (6.8), it suffices to compare the strength of the incoming 1-wave on the right of X(t) with the 
wave that corresponds on the left side: a 1-rarefaction wave if the incoming front is a rarefaction front, 
a 3-compression wave if the incoming front is a shock front. The fact that the waves on the left and 
on the right have proportional strengths is a consequence of Remark 3.1 (when the incoming wave is a 
rarefaction wave) and Proposition 3.4 (when the incoming wave is a shock). 

Finally, (6.9) is obvious since rg c Cg U C\. □ 
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Corollary 6.1. IfTV(ug) is small enough, then one has, as long as the algorithm is well- functioning, 
that 

TV(u u (t, •); (0, X(t))) + TV(u"(t, •); (X(t),L)) < CTV(u ), (6.10) 

\\u u (t, ') - Vo\\L°°(o,x(t)) + \\u u (t, •) - «o|U=°(;c(t),£) < C(TV(u ) + \\u - «olU°°(o,i))- (6-11) 

Proof of Corollary 6.1. For what concerns (6.10), we only notice that the left hand side of can be esti- 
mated by F l (t) +F 3 (t), which due to Lemma 6.2 can be estimated by CF 1 (Q). For what concerns (6.11), 
we can measure the second term in the left hand side by C(||u — Uolloo + F 1 ^)), because F 1 ^) allows 
to estimate the distance between v?(t, x), t > 0, x > X(t), and Ug(0 + ). In the same way, concerning the 
first term in the left hand side of (6.11), one has 

\u u (t,x)-u u (0+,X{0+)-)\ <CF 2 (x) fort>0, xe{0,X(t)), 

where u v (0 + , X(0 + )~) is the state on the left of the strong shock just after t = 0. Now this state 
u"(0 + ,X(0 + )-) results from the Ricmann problem (vq,Uq(0 + )). It follows from (2.28) that 

\vo-u u (0 + ,X(0+)-)\ < CK(0+)-«o| < C||«o-5o|U. 

Finally it is clear that for TV(u ) < 1, one has F x (0) < CTV(u ). D 

6.1.2 Well-functioning of the algorithm; size of rarefaction and artificial fronts 

That the algorithm is well-functioning for e small enough is a consequence of the estimates above. 

Lemma 6.3. If TV(uq) is small enough, all the states generated by the algorithm belong to T> and a 
finite number of fronts and interaction points are created. 

Proof of Lemma 6.3. If one chooses TV(uq) small enough, (6.11) ensures that as long as the algorithm 
allows to construct u", the states in u v belong to T>. Hence we know that the algorithm does not stop 
because it should generate another state. Consequently, if there is no accumulation of interaction times, 
the algorithm is functional. 

The proof that there is a finite number of fronts and interaction points resembles the case of the 
initial-value problem. New physical fronts are only generated at interaction points of weak waves with 
large amplitude (under the strong discontinuity) and at interaction points with the strong wave. But 
interaction points with large amplitude are in finite number as a consequence of (6.4). Since new 1-fronts 
under X are only generated at such interaction points, we deduce that there is a finite number of 1-fronts 
under X. Therefore there is only a finite number of interaction points with the strong wave. It follows 
that physical fronts are in finite number. 

Now we deduce that interaction points involving only physical fronts are in finite number, so artificial 
fronts are in finite number as well and finally there is a finite number of interaction points. □ 

At this stage we know that for TV(uq) small enough, the algorithm generates a front-tracking ap- 
proximation u v for all small v > 0. Let us now establish estimates on the size of the fronts that will be 
important to prove the consistency of the algorithm. 

Lemma 6.4. There exists C > such that a front a in u v satisfies: 

• if a is a rarefaction front or a compression front, then \a a \ < Cv, 

• if a is an artificial front, then \o~ a \ < Cg. 

Proof of Lemma 6.4- Consider a front a which is either a rarefaction front, a compression front or an 
artificial front. It is clear how to trace back such a front across the various interactions that it has 
undergone, to its creation locus. For that, one follows a back in time (for rarefactions and artificial 
fronts) or forward in time (for compression fronts), and at an interaction point, one follows the front 
with same nature (family and sign). In this way we trace a to its creation: 

- for a rarefaction front: at t — 0, or at an interaction point with large amplitude where no incoming 
front is of the same family, or at an interaction point on X, 
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- for a compression front: on the strong discontinuity X, 

— for an artificial front: at an interaction point with small amplitude. 

There is no ambiguity in this tracing process, since in it", rarefaction fronts of the same family do not 
interact, nor do compression fronts (which all are of the third family) or artificial fronts. 

At is creation, a rarefaction front or a compression front a has a strength \o~ a \ < is; an artificial front 
satisfies \a a \ < Cg, due to Proposition 2.2. 

Now we begin with rarefaction under X and artificial fronts. Given such a front a that one follows 
over time, one can construct 

V a (t)= Yl ( 1 + KS ik f ,S X0> x(t))Wf}\- 

f) cutting rj 
approaching a 

Reasoning as before one gets that for some C > 0, V a (i) + CQ(t) and Q(t) are non-increasing during 
the lifetime of the front a. Moreover, forward rarefaction fronts do not interact and nor do artificial 
fronts. The strength of a rarefaction front does not increase when it meets a shock in the same family 
(either the strength decreases, or the rarefaction front is killed). It follows that their strength \a a \ can 
only increase at interaction points involving a and a front of another family. At such interaction times 
V decreases and one has 

M* + )| < MO! + C'\a a (t-)\ (V a (t-) - V a (t+)). (6.12) 

It follows that t h-> |cr Q (t)| cxp (C'(V a (t) + CQ(t))) is non-increasing and that 

K{t)\ < \a a {s)\exp(C(V a (t)+CQ(t))) < \a a {s)\ exp(CTV(u )), 

which gives the conclusion for rarefaction under X and artificial fronts. 

The cases of compression fronts or rarefaction fronts above X are similar, replacing the time variable 
t with the pseudo-time variable L — x: at pseudo-times where \a a \ increases, one obtains instead of 
(6.12): 

K(x~)\ < \a a (x+)\ + C'\a a (x+)\ (V a (x+) - V a (x-)), with V a (x) ■- ^ |^|, 

/3eA a (x) 

where A a {x) :— {/3 fronts cutting C% and approaching a}. This allows to conclude as before. This ends 
the proof of Lemma 6.4. □ 

6.1.3 Total strength of artificial fronts 

Here we prove the following proposition. 

Proposition 6.1. There exists C > such that ifTV(uo) is small enough (depending only on 7, uq 
and L) and if g > is small enough (depending on v and Uq), then one has for all t > 

J2 l a «l - Cv - ( 6 - 13 ) 

a artificial front 
living at time t 

Proof of Proposition 6.1. Artificial fronts live only below the strong discontinuity, hence we only consider 
the approximation there. For this we can follow Bressan's analysis [6, 7]; we recall the argument in order 
to check that it fits in our situation. 

This analysis relies on the notion of generation of fronts. This is defined as follows: to each front in 
u v , under the strong shock, is associated a positive integer called its generation and computed by the 
following rules. Each front emerging from t = has generation 1, and when two weak fronts a and /3, of 
generation g a and gp interact: 

• if k a y^z k/j, the outgoing fronts of family k a , respectively kp, resp. k £ {k a , kp} is of generation g a , 
resp. gp, resp. max(g a ,gp) + 1, 
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• if k a — kp, the outgoing fronts of family k a , respectively k ^ k a is of generation min(g Q , gp), resp. 
m&x(g a ,gp) + 1, 

and when a 1-front a interacts with the strong 2-discontinuity, the outgoing fronts of family 3 are declared 
of generation g a . Recall that the artificial fronts are considered of family 4. 
Now we can define the functionals: 

Vfe(i) = X! 0- + ^ik a S Xa>X (t))Wa\ and Q k (t) = ^ Wa\Wfi\, 

a front cutting T t a,/3 front cutting T t 

of generation >k a approaching j3 

with max((; a ,(/^)>A; 



with k as in Lemma 6.1. Define 



V k {t) := sup Vfc(s). 
«e[o,t] 



The function Q k (t) is piecewise constant hence 5 V; consequently it can be decomposed into Q k (t) — 
Qk(t) — Qk(t) with Q k and Q& non-decreasing, and with Qfc(O) = Qfe(O). For fc > 2, one has Qfc(O) = 
and since Qk(t) > 0, one has < Qk(t) < Q k (t). 

Reasoning as in Lemma 6.1 we see that the only case when Vk can increase is when two weak fronts 
a and (3 with max(g a , gp) = fc— 1 interact, and in that case Vk{t + ) < Vu{t~) +C[Qk-i{t~) — Qk-i{t + )]- 
With Vfc(0) = for fc > 2, we deduce that for some C > 0, one has for all fc > 3, 

Vk(t)<CQ k -i(t)<CQk-i(t). 

Now Qk is modified only when a front leaves the domain (in which case it decreases) and at interaction 
times. Consider such an interaction time t involving weak fronts a and j3. Reasoning as in Lemma 6.1 
we get: 

• if mox.(g a , g ) > fc, then Q k (t + ) - Qk{t~) < 0, 

• if max( ffa , W ) = fc - 1, then Q k {t+) - Q k {t~) < Cy(r)[Q k _i(r) - Q k -i(t + )], 
. if max(.g a , g p ) < fc - 2, then Q k (t+) - Q k {t~) < CV k {t-)[Q{r) - Q{t+)}, 

When a 1-front a hits the strong 2-discontinuity, one has, due to the reflected 3-wave: 

• if 9a > fc, then Q k {t+) - Q k (t~) < C\a a \V(t~) < CV(f)[V fc (t-) - V k (t+)}, 

• if g a < k - 1, then Q k (t+) - Q k {t~) < C\a a \V k (t-) < CV k {t-)[V {t~) - V{t+)]. 
Summing these inequalities over all possible interactions making Q k increase, we get (for fc > 3): 

Q k (t) < Cii ll (0)(Q fe _i(t) + V k (t)) < CF\0)Q k ^(t). 

Hence, if TV(u ) is small enough, one has for some \i € (0, 1) that V k (t) + Q k (t) < C/J, k . In particular, 
there is some fc for which 

m < \ 

Now if initially Uq generates N fronts, then there are at most N 2 interactions of fronts of first generation, 
and hence at most C(l + -)N 2 fronts of second generation, and by induction there are at most C k (N, v) 
interaction points involving fronts of the fc-th generation at most (the point being that this function C k 
does not depend on g) . The total strength of artificial fronts of generation > fc is measured by Vj:, the 
one of artificial fronts of generation < fc is measured by C gC^^N^v) for some positive constant C. 
Hence for g small enough, the latter is less than |, which establishes (6.13). □ 
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6.1.4 Passage to the limit 

Extraction of a converging subsequence. Adding the strength of the strong discontinuity to (6.10) 
and using the definition of r, we deduce that {u v ) v> q is uniformly bounded in the space L°°(M. + ; BV(0, L)). 
Now, using that all the fronts in our approximation (including the artificial ones) have bounded finite 
speed, we classically deduce that the family is uniformly Lipschitz in time with values in L 1 (0,L): 

\\u v {t)-u v {s)\\v {Q>L) <C\t-s\ max7VK(T,-)). (6.14) 

T£[8,t] 

It follows then from Helly's compactness theorem that one can extract from (u l/ ) a converging subsequence 
(u 1 "") in L 1 locally in time and, reextracting if necessary, almost everywhere: 

u Un — >ua.e. aadinL^OjT) x (0,L)), VT > 0, (6.15) 

and the limit u belongs to L°°(IR + ; BV(0, L)) and to Cip(R + ; L\0, L)). 

The limit point u is an entropy solution. We now prove that the limit point u is a weak solution 
of the system and satisfies the entropy inequalities. For that, we first get back to conservative variables. 
We denote U v and f the functions u v and u translated in conservative variables. Using the L°° bound 
on U" 71 and Lebesgue's dominated convergence theorem, we get 

ljv„ — yjj ae andin l\{0,T) x (0,L)), VT>0. (6.16) 

Note that, since BV is an algebra, (f "") is uniformly bounded in L°°(R + ; BV(0,L)), so f belongs to 
this space as well. Using the L°°(R + x (0, L)) and Cip(R + :L 1 (0, L)) bounds on u v , we deduce that (U Vn ) 
is uniformly bounded in £ip(R + ; £ 1 (0, L)), so f also belongs to £ip(R + ;L 1 (0,L)). 

Now we consider (j],q) an entropy /entropy flux pair, with 77 convex. We include 77(f) = ±f and 
q(U) = ±/(f ) in the discussion; this will give us that f is a distributional solution. 

In order to prove the entropy inequality associated to (rj,q), it is enough to prove that for all ip £ 
£>((0,T) x (0,L)) with ip > 0, one has 

liminf J n := f [ip t ri{U" n ) + <p x q(U u ")] dtdx > 0. (6.17) 

™^+°° J(0,T)X(0,L) 

We describe the discontinuities at time t by the family of fronts {a, a € .4}; each front a has position 
x a (t) at time t, and we denote [h] a (t) the jump of the quantity h through the jump a(t). Classically we 
have, integrating by parts: 

"" dt 



J n = / y;^,^(*)){x Q (t)-[77(f "»)] a (t) -{q(U v « )]„(«)} , 

> J2 [ <P{t,x a (t)){x a (t)-[ V (U^)] a (t)-[q(U^)] a (t)}dt=: Yl J - ( 6 ' 18 ) 

a. JO a 

weak front weak front 

Here we used the fact that the strong 2-discontinuity (which travels at exact speed) satisfies the entropy 
inequality (actually, even as an equality here): 

s[v(U^)] a (t)-[q(U^)] a (t)>0. (6.19) 

This fact is general for any 2-contact discontinuity traveling at shock speed. A way to prove it is as 
follows. Denoting f + := Sk(o", f-) (with here k = 2), we differentiate 

s[rj\ - [q] = s(U_,U + ) (rj(U + ) - rj(U_)) - {q(U + ) - q(U_)) 

with respect to a and use the Rankine-Hugoniot relation to obtain 

J^( S W - [«]) = fj (v(U+)-v(U-)-Dr,(U + ) ■ (f + -f_)). (6.20) 

Here we have ff = 0, which establishes (6.19) as an equality. 

Now, let us consider the term J a in (6.18) and discuss according to the nature of the weak front a: 
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• if a is a shock, it would satisfy the entropy inequality if it was traveling at the exact shock speed; 
but since it moves at shock speed up to a small change of v n , we have in general J a > —0(l)v n \a a \. 

• if a is a rarefaction front or a compression front, one sees easily by differentiation that, s being the 
shock speed (3.9), one has 

[q(U^)} a (t)-s[ V (U^)} a (t)=0(\a a \ 2 ), 

which yields 

J a > -0(l)i/„|ff a |. 

(This could be improved in the case of compression fronts, but this has no importance.) 

Using the uniform bound on the total strength of physical waves and of artificial waves, this yields 
Jn > —Cv n , which establishes (6.17). 

The solution u is constant at some time T > 0. In our construction, after the exit time T\ (which 
satisfies (4.13)), there are only 1-rarefaction fronts in the approximation u u . Indeed, the compression 
fronts emerge from the strong 2-discontinuity only and travel backward in time. Moreover, due to (4.7) 
the rarefaction fronts travel at speed s less than 

Consider only v < ' 1 ^ '' . Hence after the time T 2 defined by 

A 2 (u ) \\\{u )\-2v 

there is no front inside the domain. Hence all the approximations u v are constant in space after some 
uniform time T 2 ; consequently so is u. 

6.1.5 Remaining cases 

We have yet to explain how we treat the cases which are not covered by (4.1). The case where (4.1) 
holds will be referred to as Case 1. The other cases are as follows. 

Case 2. Ai(u ) > and A 2 («o) > 0. This (supersonic) case is in fact by far the simplest. Indeed, in 
this case, introduce r > such that one has Ai(u) > Ai(uo)/2 on i?(uo,r). Given uq one can define on 
M the following initial data: 

uq = uq on (0, L) and uq =Tio onR\ (0, L). (6.22) 

Then if uq satisfies (1.12) with e > small enough, one can associate to this initial condition the unique 
entropy solution u in R + x R, as in [7]; moreover for e small enough, u has values in B(u ,r). For 
instance, one can obtain u as a limit of front-tracking approximations. The restriction of this solution 
to (0,Ti) x (0, L) is convenient, where 

2L 

Indeed, all the fronts have a velocity larger than Ai(uo)/2, hence leave the domain before T\, so u v (t, ■) 
is constant for all v for times t > T\. 

Case 3. A 2 (mo) < and A 3 (u ) > & Case 4. A2(wo) < and A3(uo) < 0. This cases are obtained 
from the Cases 1 or 2 above by using the change of variable x < — > —x. 

Case 5. A 2 (uo) = 0. In that case of course, Ai(uo) < and A3(uo) > 0. We let a large (and after all, 
not so large) 3-shock enter through the left side. In other words, we introduce wq — T^(a3,uo), with (73 
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negative and small. We make sure that is speed satisfies s > 3A3(uo)/4. Then we introduce the solution 
u associated to the initial data 

So = wq on R~ , uq = uq on (0, L) and uq = Uq on (L, +00). 

This solution can be constructed by the front-tracking method described above; in particular one can 
follow the 3-strong shock inside the domain by a curve X 3 (t) and get, provided that TV(uq) + \\uq — MoI|l°° 
is small enough, that the approximations satisfy 

TV(u»(t,-);(0,X 3 (t))) + TV(u v (t,-);(X s (t),L)) < CTV(u ). (6.23) 

When the 3-shock issued from x = has left the domain (0, L) (for instance at a time T = 2L/Xs(uo)), 
we are left with a state u(T, •) in (0, L) which satisfies: 

\\u(T,-)~w \\ L ~ (0 , L) +TV(u(T,-)) <KTV(u ). (6.24) 

This is proven as Corollary 6.1. Now, if |cr 3 | was small enough, one has 

Xi(w ) < and A 3 (u;o) > 0, (6.25) 

and moreover, due to 

r 3 • VA 2 = —— > 0, 

7+1 

one has 

A 2 K) > 0, (6.26) 

so we are now in position to apply Case 1. 

6.1.6 Smallness of the solution 

The last part of the proof consists in proving (1.14), provided that e > is small enough and that the 
large 2-discontinuity (and possibly the preliminary 3-shock of Case 5) is (are) well-chosen. This depends 
a bit on the cases described in Paragraph 6.1.5. 

Cases 2. &: 4. In those cases, the solution that we construct is obtained by the restriction to (0, L) of 
a solution defined on R and whose initial data has a total variation less than TV(%) + 2||u — wo||oo (see 
(6.22)). Due to Glimm's estimates, the solution satisfies that 

sup ||u(t, •) - uo\\l~(wl) + TV(u(t, •)) < K (TV(u ) + ||«o - «b|U). 
t 

Hence (1.14) follows, and here r) is a linear function of e. 

Cases 1. & 3. We only consider Case 1 by symmetry. What we have established in this case is that if 

uq =T 2 (a 2 ,Vo), a 2 < 0, 
there exists e 2 = s 2 (a 2 ) > such that for any e € (0, §2] if «o satisfies 

||«o - Bb|U-(o,L) + TV(u a ) < e, (6.27) 

then the construction given in Section 4 with this strong 2-discontinuity is valid and due to Corollary 
6.1 and the definition of r there is K = K(a 2 ) such that, including the strong discontinuity one has: 

supTV(u{t, •)) < C\a 2 \ + K{a 2 )e. (6.28) 

t 

We choose a 2 such that \a 2 \ < r//2C. Then we choose e 2 £ (0,22] such that K(a 2 )e 2 < 77/2 and we are 
done. 
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Case 5. In this case, there is a preliminary phase before getting into Case 1. In the same way as before, 
given (73 < small enough such that wo — T^ia^jUo) satisfies (6.25)-(6.26), if e is small enough and if 
Mo satisfies (6.27), then the solution that we construct satisfies 

supTV(u(t, •)) < C\a 3 \ + K'e, 

t 

||u(r, •) - w \\l~(o,l) + TV(u(T, •)) < KTV(u ), T = 2L/X 3 (u ). 

Here K' can be chosen independent of (73, by using Glimm estimates. Above, we used cancella- 
tion/correction waves for which the constant worsens as the strong shock becomes small; here this 
is not the case. 

Now, we first choose <7 3 < and Wq such that C\a^\ < rj/2. Then reasoning as before, one can find a 
size of <72 and an £2 corresponding to the second phase, with wq as a reference state, in order for (1.14) 
to be valid during this second phase. Then we find a size of £3 > corresponding to the first phase, in 
order that K' £3 < 77/2 and that the state at the beginning of the second phase is small enough to satisfy 
(1.12) with wq as a reference constant state and right hand side £2- 

6.2 Lagrangian case: proof of Theorem 4 

In this subsection, we prove Theorem 4 by adapting the arguments of Subsection 6.1 in the situation 
given by the construction of Section 5. 

6.2.1 BV estimates 

The first point is to prove uniform BV estimates on the approximation u v (again, as long as it is 
well-defined and all the states belong to V). 

For that we introduce six families of curves drawn in R + x [0, X], defined for fixed v and for fixed t 
or x. We recall that the strong 1-shock (resp. 3-shock) is represented by X\ (resp. X3), has a negative 
(resp. positive) speed, that it enters the domain at time (resp. T2) and leaves it at time T± (resp. T3). 
The curves are the following. 

• Given t £ [0, Ti], we define the curve Tj as the horizontal line segment {£} x [0, Xi(t)], glued with 
the part of curve X\ from (t,Xi(t)) to (0,L) (or a curve very close on the left of Xl). We do not 
consider Tj for t larger than T\. 

• The curve C^, defined for x £ [0,L], is obtained by gluing the vertical line segment [X^ 1 (x),T2]x{x} 
and the portion of X\ between (A{~ (x),x) and (0,L) (with in mind that this portion is "above" 
Xr). 

• The curve rf, defined for t £ [0,T 2 ], is obtained by gluing the part of the curve X\ for times in 
[t,Ti] (or a curve very close on its right), the horizontal curve {£} x [Xi(t),L] and the vertical line 
segment [t, T2] x {L}. After time T\, only the horizontal and vertical parts remain. 

• Given t £ [T2,Ts], we define the curve T^ as the part of curve X3 from (T2,0) to (t,Xs(t)) (or a 
curve very close on the right of X3) glued with the horizontal part {£} x [X 3 (t).L\. We do not 
consider Tf for t larger than T 3 . 

• The curve C|, defined for x £ [Q,L], is obtained by gluing the portion of X 3 between (T 2 ,0) 
and (X3 (x),x) (with in mind that this portion is "above" A3) and the vertical line segment 

[X^ 1 {x) 1 +^)x{x}. 

• Given t > T 2 , the curve Tf is obtained by gluing the vertical line segment [t, +00) x {0}, the 
horizontal line segment {t} x [0, A 3 (i)] and the part of the curve A 3 from (i, A 3 (i)) to (T 3 ,L) (or 
a curve very close on the left of A3). For times larger than T3, it remains only the vertical line 
segment [t, +00) x {0} and the horizontal line segment {£} x [0, L). 

We represented these six families of curves in Figure 14. 

Let us specify how we follow these curves. For i = 1,4, we just follow the curves for increasing x. For 
i = 2, we first follow the vertical line segment from top to bottom, and then the part on X\ from left to 
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(a) The curves Tj, t = 1,3,4, 6 (b) The curves Cj, i = 2, 5 

Figure 14: Curves for the By estimate in the Lagrangian case 

right. For i = 3, we first follow the part of Xi from left to right, then the horizontal line segment from 
left to right and finally the vertical line segment from bottom to top. For i = 5, we first follow the part 
of X\ from left to right and then the vertical line segment from bottom to top. Finally for i = 6, we 
first follow the vertical line segment from top to bottom, then the horizontal part from left to right and 
finally the part on X\ from left to right. In all cases we will say that we follow the curve "from left to 
right" . Given two points on one of these curves, this gives a meaning to "one is on the left of the other" . 

Now in order to get a uniform estimate on the total variation in space of u v (provided that (1.16) 
holds with e > small enough), we proceed as in Paragraph 6.1.1. For that we introduce the functionals 
for i= 1,3,4,6: 



V\t) 



E 

a front cutting T\ 



cr a \ and Q\t) 



E 



Wa\, \&p\, 



(6.29) 



a,0 front cutting VI 
a approaching 



as well as the following ones, for i = 2, 5: 

V\x)= V |<r | and Q\x) 



E 

a front cutting Ct 



E 



\<rp\ 



(6.30) 



a,/3 front cutting C^ 
a approaching /3 



Let us give some precisions for these definitions. 



• Only fronts crossing the large 1-shock on its left (resp. right) cross T^ (resp. C%, Tf) on its part 
coinciding with Xi . 

• In the same way, only fronts crossing the large 3-shock on its right (resp. left) cross Tf (resp. C®, 
r^) on its part coinciding with X3. 

• Our convention is that 2-fronts do not cut the vertical part of C% and of course that the strong 
shocks do not cut the curves. 

• If a (of family i) and (3 (of family j) cut 1^ (k — 1, 3, 4, 6), a to the left of j3, they are said to be 
approaching when i > j or i = j and at least one of a or is a shock (artificial fronts being of 
family 0.) 

• If a (of family i) and j3 (of family j) cut C% (k — 2,5), a to the left of f3, they are said to be 
approaching when: 
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— k = 2: i = l and j £ {2, 3} or i = 3 and j = 2, 

— k = 5: i = 1 and j = 3. 

Note that with respect to (6.29), we do not put a weight in the functionals V % . This is due to the fact 
that, as in this construction the strong waves are from extremal families, there is no reflected wave when 
considering the interaction "from below" between a weak front and one of the two large discontinuities. 
This allows to simplify a bit the analysis. 

Now we can get as previously the following result. 

Lemma 6.5. For C > and K > large enough, the following holds provided that TV(uq) is small 

enough. 

1. For k = 1, 3, 4, 6, the functional 

F l {t) := V l (t) + KQ l {t) is non-increasing over time. 



F (x) := V (x) + KQ (x) is non-increasing, 
F (x) := V (x) + KQ (x) is non- decreasing. 

2. One has the relations: 

F 2 (0) < CF l (T{) and F 5 (L) < C*F 4 (T 3 ), (6.31) 

F 3 (0) < F 2 (0) + F 2 (L), F 4 (T 2 ) = F 3 {T 2 ) and F 6 (T 2 ) < F 5 (0) + F 5 (L). (6.32) 

Proof. As for Lemmas 6.1 and 6.2, the first part is a consequence of Glimm's estimates for usual in- 
teractions or side interaction of weak waves; Proposition 2.2 can be applied to all those cases. Note in 
particular that at times of interaction with one of the strong shocks or at the exit time of a front, the 
functionals V l (t) and Qi(t) either decrease or stay constant (according to i and the family of the front). 
The functionals V 5 (x) and Q 5 (x) do not change when x corresponds to the location of an interaction 
with the 3-shock, and C% does not meet an exit location before x = 0. The only thing to be careful 
about concerns F 2 . Indeed there can be many "simultaneous" interaction points on C 2 when x corre- 
sponds to the position of a 2-contact discontinuity. For such a x, we analyze each interaction separately: 
making the sum of the contributions gives the same result as if the interaction times were distinct; more- 
over the 2-contact discontinuity disappears from the functionals which gives a supplementary negative 
contribution. 

For the second part, (6.31) is a consequence of Schochet's estimates for interactions with a large 
discontinuity or estimates for cancellation/correction waves at an interaction point with a large discon- 
tinuity (Propositions 3.2, 3.3, 3.5 and 3.6); moreover we notice that no front crosses (Ti,Ta) x {0} or 
(T 3 , +oo) x {L}. To get (6.32), one just has to compare the curves on which the functionals rely. □ 

One can deduce as before the following. 

Corollary 6.2. IfTV(uo) is small enough, then one has for all times t < T 2 for which the algorithm is 
well-functioning that 

TV(u v {t, ■)■ (0, X 1 (t))) + TV{u v {t, ■); (Xi(t), L)) < CTV(u ), (6.33) 

\\u v (t,-) -uo|U»(o,Xi(t)) + \\u"(t, ■) - «o"|U»(Xi(t),r,) < C(TV(u ) + \\u - uo\\l°°(o,l)), (6.34) 

where we set X\{t) = for t > T±. Moreover for all times t > T 2 for which the algorithm is well- 
functioning one has 

TV(u v (t r )\(0,X s (t)))+TV(u v (t r y,(Xs(t),L)) <CTV(u ), (6.35) 

IK(£r) -fh IU°°(o,x 3 (t)) + IK(V) ~ Vo\\L°°(x 3 (t),L) < C(TV(u ) + \\u Q -uo\\l<»(o,l)), (6.36) 
withX 3 (t) =L fort>T 3 . 

56 



6.2.2 Well-functioning of the algorithm 

Let us continue the proof by following the lines of Subsection 6.1. Lemma 6.3 is valid in the present 
situation: 

Lemma 6.6. IfTV(uo) is small enough, all the states that the algorithm generates belong to T> and only 
a finite number of fronts and interaction points are created. 

Proof of Lemma 6. 6. The first part of the statement is a direct consequence of Corollary 6.2. We focus on 
the second part. The proof has of course common features with the one of Lemma 6.3; let us nevertheless 
stress the differences. Here new physical fronts are only created at: 

• interaction points of weak waves with large amplitude, 

• interaction points involving one of the two strong shocks, 

• interaction points involving a forward 3-rarefaction R or a backward 1-compression C with a 2- 

> 
contact-discontinuity J, in the second part of the construction. This type is new with respect to 

Lemma 6.3. 

1. The proof that there is only a finite number of fronts under the strong 1-shock is actually simpler 
than in Lemma 6.3: when a weak wave interacts with the strong 1-shock, no front is reflected under 
the strong shock. Hence only interaction points of weak waves with large amplitude can increase the 
number of physical fronts under X\ ; since they make the functional F 1 (t) decrease of an amount at least 
proportional to g, they are of finite number. It follows that, under the strong 1-shock, there is only a 
finite number of physical fronts, hence a finite number of interactions with weak amplitude and a finite 
number of artificial fronts as well. 

2. Since there is a finite number of interactions with the strong 1-shock, there is a finite number of fronts 
emerging from Xi . Now between the strong 1-shock and the strong 3-shock, there is no interaction point 
that modifies the number of fronts of family 2, hence those remain finite and do not disappear before 
meeting the strong 3-shock. Let us call them from left to right (in the x variable), Ji,. . . , Jk- 

Now, define for x G [0, L] the number 

k— 1 a a 1 or 3-front, 
cutting C x , 
and on the left of Jk 

We use the same convention as before to determine when a front cuts €%. Then Af(0) is finite (since Cg 
coincides with [Ti,T 2 ] x {0} glued with X\ and no front cuts (Ti,T 2 ) x {0}). The number M{x) (for 
increasing x) can only evolve at x where C% meets an interaction point or a point where a front leaves 

the domain [0, T 2 ] x [0, L\. Through an interaction point where a front R meets a front C, N{x) actually 
stays constant. Moreover only 2-contact discontinuities leave the domain [0,T2] x [0,L] elsewhere than 
through x = L. It follows that Af(x) changes only when x corresponds to the position of a 2-contact 
discontinuity. 

Consider such a x corresponding to Jjr. Now J\f(x + ) differs from J\f(x~) for two reasons: 

• There are fronts R or C that existed before their interaction with J^ (that is for x < x), but after 
the interaction (just on the right of x) their contribution to N is 3 _fc less than before. 

• Each interaction point on J^ generates a new front of type R (resp. C) when the incoming front 
(the corresponding front at x < x) is of type C (resp. R). The contribution to J\f(x + ) of such a 
new front is T,k=k+i 3 _fc < 2 • 3"*" 1 . 

It follows that J\f is non-increasing in [0, L] and Af loses at least 3~ K through each interaction point 

involving a J front that it meets. Hence the number of interactions between a front R or C and a J 
front is finite and consequently so is the number of fronts between the two strong shocks. 
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3. Since there is a finite number of fronts above the strong 1-shock, these generate a finite number of 
interactions with the strong 3-shock. Consequently a finite number of fronts emerge from the strong 
3-shock. But above this strong 3-shock, there is no longer any creation of fronts. In final, there is a finite 
number of fronts, and consequently of interaction points. □ 

6.2.3 Conclusion 

The rest of the proof is very close to Subsection 6.1. In particular, the statement of Lemma 6.4 on the 
size of rarefaction, compression or artificial fronts is valid as it stands. The proof can be adapted without 
difficulty and consequently we omit it. The same is true for Proposition 6.1 regarding the total strength 
of artificial fronts and the same argument can be used (even, a bit simplified by the absence of reflected 
waves); again there is no artificial front above the strong shock. Next the arguments allowing to pass to 
the limit and obtain a weak entropy solution u can be entirely reproduced from Subsection 6.1 except for 
the proof that the strong shocks satisfy (6.19). Here we use that in (6.20) the second factor is positive 
(by convexity of n) and that the shock speed is increasing along the shock curve (see (3.25)). 
It is finally sufficient in order to conclude to prove the following. 

Lemma 6.7. For a time T4 satisfying 

n<n+ 2 ^ (6.37) 

|Ai(« x )|-2i/ 

one has for all v > small that 

u u (t, •) is constant for t > X4. 

Proof of Lemma 6. 7. Above the strong 3-shock there are no 2-contact discontinuity, but only fronts of 

type R and C . Since no new fronts are created, backward compression fronts C live only before time T3. 

It remains to consider the fronts R, which emerge before time T3, and travel through the domain from 
right to left at speed at least of — (Ai(uf )/2) — v. The conclusion follows. □ 

The last part of the proof of Theorem 4 consists in proving (1.18). Again, we can prove that that 
given Vq and wj~ satisfying the requirements of Subsection 5.1, there exists £0 such that for any e < e , if 
uo satisfies (1.16), then the construction above is valid. Since Vq and v± can be chosen arbitrarily close 
to Ho, the conclusion follows as in Paragraph 6.1.6. 

6.3 Moving to any constant 

6.3.1 Eulerian case: end of the proof of Theorem 1 

As we mentioned earlier, once the state of system is driven to a constant, reaching any constant in the 
Eulerian situation (keeping an arbitrarily small total variation in x uniformly in i) can be seen as a 
consequence of [32] and [16]. Indeed, a corollary of these results is the following. 

Theorem 5 ([32, 16]). Consider system (1.1). For any u* £ D, and any n > 0, there exists 6 > and a 
time T such that, for any u a ,Ub G C rl ([0, L]) satisfying 

\\u a - u*|| c i([ ,l]) < $ and \\u b - u*||ci([o,z,]) < 8, 
there exists a C 1 solution u of the system driving the state from u a to Ub, with |Jm||c°([o,t];C 1 ([o,l])) < V- 

Actually [32] allows to treat the case where ^ {Ai(u*), \ 2 (u*), A3 (it*)} and [16] the remaining cases, 
relying on r\ ■ VA2 7^ at points where A2 = and on r 2 ■ VA; 7^ at points where A, = 0, i = 1, 3. 

This statement applies of course to u a and Ub constant, but it is local. Now we deduce a global result: 
given ui,u~2 € fi, let us explain how to drive u\ to u 2 - We consider a smooth curve j(s) from ui to U2- 
For 77 > 0, the above statement gives us a S(s) for each point u* := 'f(s) of this curve. By compactness 
of the curve, we can extract a finite (sub)cover of 7 by balls B(7(sfc), <5(sfc)/2). Then one can drive from 
u\ to u~2 by successive steps leading a j(sk) to another. The resulting solution has constantly a C^norm 
in space less than n. This ends the proof of Proposition 1.1 and hence of Theorem 1. 
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6.3.2 Lagrangian case: end of the proof of Theorem 2 

We begin by proving Proposition 1.2. Here we cannot apply [16] to treat the vanishing characteristic 
speed A2, because it is identically zero. 

The first step of the proof is to go from the constant state u a to another constant state u' for which 
S(u') — S(ub). This relies on the following. 

Lemma 6.8. Let Uo = (to, vq, Pq) £ 51 and 77 > 0. For any x > S(uo), there exists T > and an entropy 
solution u with initial data Uo, such that 

Vte[0,T], TV{u{t, •)) < T) and S(u{T)) = X - (6-38) 

Proof of Lemma 6.8. Of course this could not be achieved via regular solutions. Here the idea is to 
use a succession of 1-shocks and 1-rarefactions. We use the parameterizations of Si and IZi by x as in 
Paragraph 3.1.2. 

Starting from some u = (t^v^P), one introduces u :— S\(x,-)u and u := 1Z\{\/x,-)u. One can 
construct an entropy solution from u to u by letting first the shock (u, it) cross the domain from right to 
left during some interval [0, T{\ and then letting the rarefaction wave (u, u) cross the domain from right 
to left during some interval [Ti,T 2 ]; both waves have of course a negative speed. We write u — (f, w, P) 
and u — (t,v,P). Using formulas (3.22) and (3.26), we obtain that: 



P = xP and P 




It follows that 

TV{u{t, ■)) <P(x-l)+T (j^j - 1 ) + VPt^ - " i ^_^ [or/^lo.r,.]. (<,.:■!<)) 



TYinU.-U :''(.r-l)IT^(.r l 1) + _ ^ ^lx^-^1 - x~<). for i € [Ti , T 2 ] . (6.40) 

It is clear that the right hand sides go to zero as a; — > 1 + . Consequently there exists x £ (1,2] depending 
on 77 such that, for any x £ [1,3?], whenever 

( e x \ 1/7 - 

r< ((7-1)^) and P < P a , (6.41) 

the corresponding solution u satisfies 

TV{u{t,-))<n. (6.42) 

We notice that 

P + x^ 1 



q(x) :— x — is increasing as a function of x, 

\j3x + i; 



so that g(x) > 1 for x > 1. We introduce n by 

X ~ S(u ) 

.10g(£/(x)) 



(6.43) 



Now we obtain a solution to the claim by letting successively n shocks (u 1 , u l ) /rarefaction waves (u l , it, 1 ) 
cross the domain, with u° := Uo, u l — Si(x,u l ), u % — TZi(l/x,u l ) and u l+1 = u l , i — Q,...,n — 1; 
and then on lets a last shock (u n , u n+1 ) and a last rarefaction wave (u n+1 , u n+1 ) cross the domain with 
u n+1 = 5i(a;',u n ) and u n+1 = Tli(l/x', u n+1 ) where x' is such that 

log(s(a/)) = x - S'(uo) - nlog(g(x)). 
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Note that clearly x' £ [1, a?] and that after the passage of this last shock, the state has reached the 
entropy %. 

Let us finally justify that the states appearing during the passage of the successive shocks/rarefaction 
waves satisfy each (6.42). It is easy to see that at each stage the intermediate state u = (Tk,Vk,Pk) 
satisfies 

T k = T g(x) k ^, P k = P and S(u k ) = S(u ) + klog(g(x)) < X - 



We deduce that 



pS(u k ) x 

(7-l)^-<(7-l)p-. 

Pk r o 



Hence the state u satisfies (6.41). It follows that (6.38) is satisfied. □ 

Remark 6.1. The rarefactions above do not change the physical entropy S; however they are useful to 
ensure that the pressure does not become large, which would be costly in terms of total variation of the 
solution u. 

End of the proof of Proposition 1.2. We apply Lemma 6.8 to drive the state u a to a state v! = (r', v',P') 
for which S(u') — S(u b ). It remains to drive u' to u b — (r b , v b , P&). Now, for the isentropic Euler equation 
in Lagrangian coordinates: 

dtT — d x v = 0, 

d t v + d x P = 0, (6.44) 

P = S{u b )T-\ 

one can find a time T > and a C 1 solution (r, v) satisfying 

(t,iO|*=o = {t',v') and (T,v)\ t=T = (n,V b ), 

and 

Vt£[0,T], TV((T,u)(t r ))<n. 

Indeed, since here the characteristic speeds do not vanish at all, it is a consequence of [32] that Theorem 5 
is valid for System (6.44) (see also Gugat-Leugering [26] for a related result). As before Theorem 5 gives 
a local solution, but one can reason as in Paragraph 6.3.1 to drive the solution along a curve from 
(t',v') to (rfe,Vb). Now, this regular solution (t,v) of the isentropic model gives a fortiori a solution of 
the non-isentropic model by setting P — S'(ub)r^ 7 . And this solution drives v! to u b as required, with 
C°([0, T]; C 1 ([0, L])) norm bounded by Cr\. Reducing r\ if necessary, this gives a solution to the problem; 
this ends the proof of Proposition 1.2. □ 

Remark 6.2. The fact that solutions of (6.44) give particular solutions for (1.3) is true for regular 
solutions but fails for weak solutions. See Saint-Raymond [37] for more details. 

End of the proof of Theorem 2. Consider TLq and u\ as in the assumptions of Theorem 2. We apply 
Theorem 4 to TLq and with S > sufficiently small to ensure that 

i - | ^ x CM . S(ui) + S(u ) 
\u - u Q \ < 6 =*> S{u) < . 

We find an e > 0, a T > and a solution driving uo to some constant state u\/2 £ B(uq, S). Now with 
Proposition 1.2 we drive U1/2 to u\ with a sufficiently small C°([0, T]; C 1 ([0, L])) norm, and this gives a 
suitable solution to the problem. □ 

7 Two final remarks 

7.1 About the dependence of e with respect to rj 

Here we give an idea of the size of e with respect to rj. Only the first phase of our construction (driving 
uq to some constant) is of interest, since for the second phase, the total variation in space can be chosen 
arbitrarily small uniformly in time, independently of e. 
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Eulerian case. We first consider the Eulerian case. Discussing as in Paragraph 6.1.6, we see that it is 
enough to estimate K{a2) as a function of oi (Cases 1, 3 and 5). Recall that the constant K' appearing 
in the Case 5 is independent of the strength a 3. As we explained, the other two cases (2 and 4) give an 
e depending linearly of r\. 

Now the coefficients of interaction with a strong shock corresponding to K(7j2) are of two types: 
either coefficients corresponding to standard interaction coefficients (these are bounded as <7 2 tends to 
zero), and those who correspond to the use of Proposition 3.4. The latter coefficients are not bounded 
as <72 tends to zero, but it is not difficult to check from (3.46) that K(a2) is of order l/l^l- To make 
the right hand side of (6.28) less than 77, one takes 02 of order 77 and then e such that K(a 2 )£ < f}/2. 
This involves that e is of order y/rj. 

Lagrangian case. Now, for what concerns the Lagrangian case, the coefficients corresponding to the 
use of the two large shocks are of order K(ai) = (D(l/\(fi\ 2 ), i = 1, 3, sec Propositions 3.5 and 3.6. This is 
due to the fact that in this case we use interactions within a family to get correction/cancellation waves. 
An important fact is that the fronts that are cancelled along the second shock (that is, the 3-shock), have 

a total variation or order 0(l)e, not 0(e/\ai\ 2 ). Indeed, they are all J waves which were generated either 

by a simple interaction with the strong 1-shock (in the case referred to as Group I in Subsection 5.2) or 

> 
by Proposition 3.5 (in the case referred to as Group II). But in Proposition 3.5, the outgoing J wave 

has the same order of strength as the incoming weak wave independently of ~o\ (see in particular (3.51), 

(3.52), (3.53) and Remarks 3.5 and 3.7) - this not the case for the strength 71 of the additional C fan. 
This involves that here e is of order Nf\. 

7.2 About the time of controllability 

We conclude with an informal discussion about the time of controllability, in particular for small ?/. It 
is clear that when 77 is small, it is costly in terms of time of controllability: think for instance about the 
case where ||zZi — wolloo 3> 77. 

Lagrangian case. We begin with the case of System (1.3). In that case, the time of the first phase 
(driving uq to some constant) is bounded easily using (5.12), (5.13), (5.14) and (6.37) by: 

T<2L ( l 2 



|Ai(«b)l AaOtf) |Ai(«r)| 

As 77 tends to zero, one can in fact let the "2L" be closer to L (see (5.6), (5.7), (5.8) and (5.9)), and the 
characteristic speeds above converge to their values at uq. It follows that as 77 — > + , one can estimate 
the time of this first phase as 



|Ai(u )| A 3 (u ), 

In particular, the time of controllability in this first phase is not affected by 77. 

Concerning the second phase, that is, driving the constant state U1/2 to the final constant state u\, 
there are two parts. The second one consists in using a C 1 solution to drive from v! to Ui by using a 
C 1 solution of the isentropic equation. It is not difficult to see that one needs 0(\u\ — u'\/rj) steps for 
this (for instance, one uses rarefaction waves and regular compression waves). The first part is more 
expensive. Indeed, to go from M1/2 to vf, one uses n steps, with n defined in (6.43). But one can see that 

9(*) = * f# ± ^V = 1 + i^i* I) 3 + 0((x l) 4 ), 



f3x + 1 J 12 7 3 

while (6.39)-(6.40) gives a total variation of order x—1. It follows that here n = 0((S(ui) — S(ui/2))/r] 3 ). 
Moreover, as r\ — > + , u\/2 gets closer to uq, so one can roughly estimate the cost of this second phase as 

77 d 
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Eulerian case. In the case covered by (4.1) (or its vertically symmetric) the time of controllability to 
some constant (that is, of the first phase) is estimated by (6.21). Then we let v go to 0, and as r\ — > + , 
the 2L gets closer to L (see again the definition of r in this case); hence one can estimate the time of 
controllability in this first phase by 

T ~ 



A2(«o) |Ai(u )| 

In the symmetric case (Case 3), one replaces Ai with A3 and puts absolute values on A2. In the supersonic 
Case 2 (resp. Case 4), it is easy to check that 

L 

resp. 



Ai(uo) V ' |A3(wo)| 

The critical Case 5 is more complex. One can use a 3- wave to shift the characteristic speed A2, but as 
this wave is of order rj, the resulting A2 is of order r\ as well; it follows that in this case the time of the 
first phase depends on rj and can be estimated in the rough form 

T = o(- 
\V 

Now, concerning the second phase from U1/2 to u\ 1 reasoning as before, one can see that, if there is no 
critical state (making a characteristic speed vanish) on the way, this second phase needs C(|wi — W1/2I/77) 
steps, so 

Tl = o nu 1 -u \ \ (7l) 

When at some point a characteristic speed A^ vanishes, we use waves of the other families to drive the 
state of the system away from the characteristic manifold {A^ = 0}. This costs a time of order 0{k/tj) 
to obtain |Ai| > k. Then one can move in the direction r^ by letting successively i-rarefaction waves or 
i-compression waves cross the domain. Each step in the direction r.; costs \j k in time for a displacement 
of order r\. Hence to obtain a displacement of ari in this context, we use a time of 

\Kri rj 

which indicates that it is favorable to take k = -y/a, and this gives a cost of order O ( > Ul ~ u °> — j . 
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